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Abstract 

This paper describes a topological search for an invisibly decaying Higgs boson, H, produced via 
the Bjorken process (e + e~ — ► HZ). The analysis is based on data recorded using the OPAL de- 
tector at LEP at centre-of-mass energies from 183 to 209 GeV corresponding to a total integrated 
luminosity of 629 pb -1 . In the analysis only hadronic decays of the Z boson are considered. A 
scan over Higgs boson masses from 1 to 120 GeV and decay widths from 1 to 3000 GeV revealed 
no indication for a signal in the data. From a likelihood ratio of expected signal and Standard 
Model background we determine upper limits on cross-section times branching ratio to an invis- 
ible final state. For moderate Higgs boson decay widths, these range from about 0.07 pb (Mg = 
60 GeV) to 0.57 pb (Mh = 114 GeV). For decay widths above 200 GeV the upper limits are of the 
order of 0.15 pb. The results can be interpreted in general scenarios predicting a large invisible 
decay width of the Higgs boson. As an example we interpret the results in the so-called stealthy 
Higgs scenario. The limits from this analysis exclude a large part of the parameter range of this 
scenario experimentally accessible at LEP 2. 
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1 Introduction 



An intense search for the Higgs boson was undertaken by all of the four LEP experiments in 
various Standard Model and non-Standard Model search channels. Searches for the Standard 
Model Higgs boson, exploiting the prediction for its decay modes and also searches for invisible 
Higgs boson decays as predicted by various extensions of the Standard Model excluded Higgs 
masses up to 114.4 GeV [1, 2]. These latter searches assumed a rather small invisible decay width 
comparable to the predicted Standard Model decay width for a light Higgs boson and well below 
the experimentally achievable mass resolution of about 3 to 5 GeV. 

Recent theories that postulate the existence of additional spatial dimensions offer a new 
possibility for invisible Higgs decays [3]. In such theories the Planck mass is lowered to the TeV 
range and a rich spectrum of new particles appears, like graviscalars in the case of flat extra 
dimensions. Hence the Higgs boson can mix with the graviscalars, which leads to a missing 
energy signature in the detector [3]. This mixing can result in a large invisible decay width of 
the Higgs boson, depending on the model parameters, and would therefore alter the Standard 
Model branching ratios. As a consequence of the broadening of the Higgs resonance in the 
recoil mass spectrum, the signal-to-background ratio can deteriorate significantly. In a worst 
case scenario, searches optimised under the assumption of a narrow Higgs resonance might have 
missed the detection of a kinematically accessible Higgs boson at LEP. 

This paper describes a search for the Higgs boson, H, which imposes no constraints on the 
total decay width. The Higgs boson is assumed to be produced in association with a Z boson 
via the Bjorken process, e + e~ — > HZ, where the Z is required to decay hadronically and the 
invisible Higgs boson is detected as missing energy Emis in the event. The results are presented 
in a model-independent way in terms of limits on the Bjorken production cross-section times 
branching ratio, a(Mn,Tn) x BR(H — > Emis), at a centre-of-mass energy of 206 GeV, where 
Th is the Breit-Wigner width of the Higgs boson. A simple model extending the Standard 
Model with additional SU(3)cxSU(2)lxU(1)y singlet fields which interact strongly with the 
Higgs boson ("stealthy Higgs scenario" [4]) is chosen as an example for the interpretation of the 
result. This interaction gives rise to a large invisible decay width of the Higgs boson. This 
dedicated search expands on the previous decay-mode-independent search [5] carried out by the 
OPAL Collaboration which reported for the first time limits on the HZ production cross-section, 
interpreted in the stealthy Higgs model. 

The paper is organised as follows. Section 2 introduces the stealthy Higgs scenario. Section 
3 gives details about the modelling of signal and background. Section 4 describes the search and 
the results are interpreted in Section 5. We summarise the results in Section 6. 

2 The stealthy Higgs scenario 

In general renormalisable theories there might be other fundamental scalars, in addition to the 
Standard Model Higgs boson, that do not interact with normal matter. To investigate the influ- 
ence of a hidden scalar sector on the Higgs observables the stealthy Higgs scenario conjectures 
the existence of additional SU(3)c><SU(2)lxU(1)y singlet fields called phions. Radiative cor- 
rections to weak processes are not sensitive to the presence of singlets in the theory because no 
Feynman graphs containing singlets appear at the one-loop level. Since effects at the two-loop 
level are below the experimental precision, the presence of a singlet sector is not ruled out by any 
LEP 1 precision data [4]. These phions would not interact via the strong or electro- weak forces, 
but couple only to the Higgs boson [4], thus offering invisible decay modes to the Higgs. The 
width of the Higgs resonance can become large if either the number of such singlets, N, or the 
coupling, lo, is large, thus leading to a broad mass spectrum recoiling against the reconstructed 
Z boson. 
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Neglecting the fermionic part for the moment, the Lagrangian of the scalar sector in this 
model, adding only four new parameters to the Standard Model, is given by 



•^scalar — -^-Higgs ~\~ -^-phion -^interaction (1) 

£ffi g gs = -d^^-M^-y) 2 (2) 

11 K 

£phion = -2 9 l*<P &i <P ~ 2 m phion ^ _ giV ^ ^ ^ 

-interaction = ~ I %= 2 $</>■ (4) 



2VN 

The term -Cffiggs describes the usual Standard Model Higgs doublet 4> acquiring the Standard 
Model vacuum expectation value, v, and having its self-coupling A. In the free Lagrangian of 
scalar singlets, £ p hiorn the singlets with mass m p hi n are denoted as the 0(-/V)-symmetric multiplet 
(p. The phions also have a self-coupling k, which is fixed at k{2Mz) = 0, to allow for the widest 
parameter range of the model. The self coupling term entering loop calculations is suppressed like 
1/N. The interaction term between the Higgs and the additional phions, interaction, leads to the 
phenomenological consequence of invisible Higgs decays because the Higgs boson couples to the 
phions independently of their mass. The strength of the coupling is instead proportional to the 
coupling constant u, which is a free parameter of the model. Even though the vacuum-induced 
mass term of the phions after the symmetry-breaking is suppressed like 1/y/N [4], the phions 
occur in loop corrections to the Higgs boson propagator and therefore affect the resonance width 
of the Higgs boson. An analytic expression [6] for the change in the Higgs width compared to 
the Standard Model decay width, Tsm, can be found in the limit N — > oo, when neglecting the 
self-coupling of the phions as a small effect: 

,2„,2 



r H (M H ) = r SM (M H ) + x ^1 - Am^JMl (5) 

The cross-section for the Bjorken process can be calculated from Equations 9 and 10 of refer- 
ence [4]. Using the parametrisation for the invisible decay width (Equation 1 and 2 in [7]) one 
can express the total cross-section for the production and invisible decay by 



/^J~gJ pinv 
dsi a {e+e -^zn)(si) ff((M 2_' a/)2 + fl/r 2 r (6) 

Here si denotes the invariant mass squared of the invisible decay products of the Higgs boson. 
The production rate of these invisible masses is given by the Standard Model cross-section 1 
C( e +e--»ZH)( s -0 f° r a Higgs boson of mass y/sj. Hence the Standard Model cross-section com- 
pletely determines the dependence of the total cross-section on the centre-of-mass energy (see e.g. 
in [8]). Therefore the total cross-section goes rapidly to zero for Higgs boson masses above the 
kinematic limit. The effect of the convolution with the Breit- Wigner-like function is a broadening 
of the resonance in the recoil mass spectrum and hence a dilution of the signal-to-background 
ratio. In extreme cases of large invisible decay width one could expect the Higgs recoil mass spec- 
trum to mimic the background. In such extreme cases even a light and kinematically accessible 
Higgs boson might have escaped detection at LEP. 

In Section 5.2 we derive limits on the stealthy Higgs model under the assumption of m p hion = 
0. By simulating signal spectra for different Higgs boson masses Mh and widths Th, we set limits 
in the cj-Mh plane in the large N limit. 

x By choosing uj > one can prevent the phions from acquiring a non-zero vacuum expectation value and avoid 
a Higgs-phion mixing due to a non-diagonal mass matrix. In case of non-zero mixing, the couplings of the lightest 
scalar to the gauge boson would decrease proportional to the cosine of the mixing angle. As a consequence the 
cross-section of the Bjorken process would be lowered. 
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3 Data sets and Monte Carlo samples 



3.1 The OPAL detector and event reconstruction 

The OPAL detector [9], operated between 1989 and 2000 at LEP, had nearly complete solid 
angle 3 coverage and excellent hermeticity. The innermost detector of the central tracking was a 
high-resolution silicon microstrip vertex detector [10] which lay immediately outside the beam 
pipe. The silicon microvertex detector was surrounded by a high precision vertex drift chamber, 
a large volume jet chamber, and z-chambers which measured the z coordinates of tracks, all in a 
uniform 0.435 T axial magnetic field. A lead-glass electromagnetic calorimeter with presampler 
was located outside the magnet coil. In combination with the forward calorimeters, a forward 
ring of lead-scintillator modules (the "gamma catcher"), a forward scintillating tile counter [9, 11], 
and the silicon-tungsten luminometer [12], the calorimeters provided a geometrical acceptance 
down to 25 mrad from the beam direction. The silicon-tungsten luminometer served to measure 
the integrated luminosity using small angle Bhabha scattering events [13]. The magnet return 
yoke was instrumented with streamer tubes and thin gap chambers for hadron calorimetry and 
is surrounded by several layers of muon chambers. 

The analysis is based on data collected with the OPAL detector at LEP 2 from 1997 to 
2000 at centre-of-mass energies between 183 and 209 GeV. The integrated luminosity analysed 
is 629.1 pb" 1 . To compare with the Standard Model Monte Carlo the data are binned in five 
nominal centre-of-mass-energy points, corresponding to the energies at which the Monte Carlo 
is produced, as detailed in Table 1. 

A fast online filtering algorithm classifies the events as multi-hadronic. Events are recon- 
structed from tracks and energy deposits ("clusters") in the electromagnetic and hadronic calorime- 
ters. All tracks and energy clusters satisfying quality requests similar to those described in [14] 
are associated to form "energy flow objects". The measured energies are corrected for double 
counting of energy in the tracking chambers and calorimeters by the algorithm described in [14]. 
Global event variables, such as transverse momentum and visible mass, are then reconstructed 
from these objects and all events are forced into a two-jet topology using the Durham algorithm 
[15]- 

3.2 Signal and background modelling 

To determine the detection efficiency for a signal from an invisibly decaying Higgs boson and the 
amount of expected background from Standard Model processes, several Monte Carlo samples 
are used. Signal events for a hypothetical Higgs boson mass Mh decaying with arbitrary broad 
width Th are simulated by reweighting invisibly decaying events of type H — ► XiXi- The mass 
of neutralinos x? 1S chosen such that the Higgs boson with mass rrii can decay into a pair of 
neutralinos, which leave the detector without being detected. These Higgs bosons with decays 
into 'invisible' particles are generated with masses rrii from 1 GeV to 120 GeV with the HZHA [16] 
generator. The HZHA events are generated assuming the Standard Model production cross- 
section <7( e + e -^ ZH ) for the Higgs boson. The test masses rrii are spaced in steps of 1 GeV. The 
spacing of the test masses is chosen such that they are not resolved by the detector in the signal 
yielded after a reweighting procedure described in the following. From Equation 6 one extracts 
the event weights Wi{rrii; Mh,Th) for a mass point rrii contributing to the search for a Higgs 
boson of mass Mh and total decay width Th- The total decay width Th is defined as the sum of 

3 OPAL used a right-handed coordinate system. The z axis pointed along the direction of the electron beam 
and the x axis was horizontal pointing towards the centre of the LEP ring. The polar angle 9 was measured with 
respect to the z axis, the azimuthal angle <j> with respect to the x axis. 
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Standard Model width and invisible width . 



Wi(mf, Mh, Th) = 120Ge v l^T—, ( ? ) 

2^m J= lGeV dnij \ m 3> 



drrii 



da ( ff(e+e-^ZH)(m t ) 2mfrg v 

dm^ 7r((M 2 -m 2 ) 2 + m 2 r 2 )" l8j 

The Standard Model cross-section 0"( C + C -^ ZH ) for the Bjorken production process in Equation 8 
propagates the centre-of-mass energy dependence of the total cross-section into the weights. The 
unweighted signal Monte Carlo samples contain 2000 events per mass point m^. In the reweighted 
signal Monte Carlo sample all test masses contribute according to their weight. The reweighted 
masses Mh range from 1 to 120 GeV spaced in steps of 1 GeV. The smallest width simulated 
by this procedure is a Th of 1 GeV and the largest a Th of 3TeV. The detection efficiency for 
a Higgs boson with Mh and Th is estimated by the sum of selected event weights assuming 
binomial errors. 

The classes of Standard Model background processes considered are two-photon 2 , two- and 
four-fermion processes. For simulation of background processes the following generators are 
used: KK2F [17] and PYTHIA [18] (qq(7)), GRC4F [19] (four-fermion processes), PHOJET [20], 
HERWIG [21], Vermaseren [22] (hadronic and leptonic two-photon processes). For Monte Carlo 
generators other than HERWIG, the hadronisation is done using JETSET 7.4 [18]. The integrated 
luminosity of the main background Monte Carlo samples is at least 15 times the statistics of the 
data for the two-fermion background, 24 times for the four-fermion background and 30 times for 
the two-photon background. The Monte Carlo events are passed through a detailed simulation 
of the OPAL detector [23] and are reconstructed using the same algorithms as for the real data. 



4 Search for e + e — > HZ with Z — > qq and H — > .Emis final state 

The event selection is intended to be efficient for the complete range of possible Higgs masses 
Mh and corresponding decay widths Th studied in this search. The preselection cuts remain 
relatively loose because the final discrimination between signal and background is done by a 
likelihood-based selection. The optimised likelihood selection has to account for the fact that 
the kinematical properties of the signal change considerably over the range of masses and width 
hypotheses considered. 



4.1 Event topologies 

The signal signature is generally characterised by an acoplanar two-jet system from the Z boson 
decay. We use the term 'acoplanar' for jet pairs if the two jet axes and the beam axis are 
not consistent with lying in a single plane. The decay products of the Z boson are preferentially 
emitted into the central part of the detector, recoiling against the invisibly decaying Higgs boson. 
This is because, in contrast to the irreducible background of ZZ — > qq vv which is produced 
with an angular dependence of the differential cross-section proportional to cos 2 9, the Bjorken 
process is proportional to sin 2 6. The Higgs boson decay leads to a large missing momentum and 
a significant amount of missing energy. In two-photon processes, where the incoming electron 
and positron are scattered at low angles, usually one or both of the electrons remain undetected. 
Events of this type have large missing momentum with the missing momentum vector, puis, 

2 Two-photon interactions occur when an electron and a positron at high energies and in close proximity emit 
a pair of photons which interact via the electromagnetic force to generate a fermion pair. 
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pointing at low angles to the beam axis. The two-photon events have a small visible invariant 
mass Myis and a tiny transverse momentum p^is but a considerable longitudinal momentum 
along the z-axis in the common case that the two photons do not have equal energy. Due to 
these special characteristics this background can be easily reduced to a negligible level. 

The two-fermion background important for this search consists of Z/7* — > qq(7) events. 
These events tend to have a big cross-section if one or more initial state radiation photons 
(abbreviated as ISR photons) are emitted so that the effective centre-of-mass energy is 
reduced to a value near the Z-resonance (so-called radiative return events). The emission of ISR 
photons happens predominantly at small polar angles. In case of a mismeasurement or escape 
of the ISR photons through the beam pipe these events have a sizeable missing momentum 
preferentially oriented at small polar angles, close to the beam pipe. In such events the two jets 
are almost coplanar. 

The most difficult background to separate is four-fermion processes with neutrinos in the 
final state, such as W + W - ^ l^v qq and W^e^v — ► qqe^f with the charged lepton escaping 
detection. The irreducible background to this search stems from ZZ — ► z^qq (about 28 % of all 
ZZ decays) leading to a signature indistinguishable from a signal event with a Higgs mass close 
to the Z boson mass. The vector bosons are usually not produced at rest, leading to a transverse 
momentum of the two-jet system and therefore to a large acoplanarity of the jets, as in the signal 
case. Furthermore, the missing momentum vector points into the central detector more often 
than for the two-fermion case. To discriminate between this background and the signal one can 
exploit the difference in the angular distribution of the differential production cross-section. 

4.2 Preselection 

The final discrimination between signal and background is made using a likelihood selection. 
Therefore the preselection cuts remain relatively loose. From cut No. (5) onwards the cut values 
have been tuned at a centre-of-mass energy of 206 GeV using a simple figure of merit (FoM) that 
is independent of the a priori unknown signal cross-section. Given only the signal efficiency e 
and the background expectation B exp as a function of a cut, this figure is defined as FoM = 
2.1 x \/B cx p + 2.4/e. This FoM gives a reasonable approximation to the median number of signal 
events that can be excluded at 95 % confidence level using the likelihood ratio. Hence a smaller 
FoM corresponds to a better preselection. 

In order to reduce the amount of background data only a subset fulfilling the following quality 
criteria is analysed. The following cuts remove almost all the two-photon background: 

(1) To reduce two-photon and accelerator induced background, track criteria are applied de- 
manding that more than 20 % of all tracks be qualified as good measured tracks [24] and 
that at least 6 of them be found. 

(2) A forward energy veto rejects events with more than 5 GeV in either the left or right 
compartment of the gamma catcher calorimeters or the silicon tungsten luminometers. 
Events with more than 2 GeV in the forward calorimeters are also removed. 

(3) The missing transverse momentum p^jis should exceed 1 GeV and iWyiS has to be larger 
than 4 GeV. 

(4) Less than 20 % of the measured visible energy .EviS should be located close to the beam 
pipe in the region | cos#| > 0.9. 

(5) The visible energy -Eyis must be less than 90 % of \/s. 
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(6) It is required that the visible mass of the event should be of order M%, i.e. 55GeV < 
Myis < 105 GeV. An asymmetric cut around the Z mass is chosen, since with increasing 
Higgs mass Mr the Z bosons will be more and more off-shell. 

The remaining backgrounds at this stage, which are more difficult to remove, are mismeasured 
Z/7* — > qq events, four-fermion processes with neutrinos in the final state, such as W + W" — > 
i^vqq and W^e^u — > qqe T u with the charged lepton escaping detection (see Table 2). 

(7) To select events that are well measured in the detector with a visible mass Myis close to 
Mz and a sizeable transverse momentum py Ig the following criterion is applied: Myis + 
5 x p^ ls > y/s/2 

(8) A large part of the qq events and the remaining two-photon background is eliminated by 
requiring the visible transverse momentum j>y IS > 6 GeV. 

(9) To remove backgrounds in which particles go undetected down the beam pipe, the pro- 
jection of the visible momentum along the beam axis, p V iS> 1S required to be less than 
0.294^. 

(10) To reduce the radiative qq(7) background, the polar angle of the missing momentum vector 
must lie within the region | cos#mis| < 0.9. 

(11) The axes of both jets, reconstructed with the Durham algorithm, are required to have 
a polar angle satisfying |cos#| < 0.9 to ensure good containment. Furthermore this cut 
exploits the fact that events of the WW and ZZ background are produced according to an 
angular distribution proportional to cos 2 #. 

(12) The remaining background from Z/7* — > qq is characterised by two jets that tend to be 
back-to-back with small acoplanarity angles, in contrast to signal events in which the jets 
are expected to have some acoplanarity angle due to the recoiling Higgs boson. Here the 
acoplanarity angle <^>acoplan is defined as 180° — <pjj where 4>jj is the angle between the 
two jets in the plane perpendicular to the beam axis. This background is suppressed by 
requiring that the jet-jet acoplanarity angle be larger than 5°. 

(13) W + W~ events with one of the W bosons decaying leptonically and the other decaying 
into hadronic jets are rejected by requiring that the events have no isolated leptons. In 
this context, leptons are low-multiplicity jets with one, two or three tracks, associated to 
electromagnetic or hadronic energy clusters, having an invariant mass of less than 2.5 GeV 
and momentum in excess of 5 GeV. In the case of a single-track candidate, the lepton is 
considered isolated if there are no additional tracks within an isolation cone of 25° half- 
angle, and if the electromagnetic energy contained between cones of 5° and 25° half-angle 
around the track does not exceed 5 % of the sum of the track energy and the electromagnetic 
energy within the 5° half angle cone. In the case of a two- or three-track candidate, 
consisting of the tracks and electromagnetic or hadronic energy clusters confined to a cone 
of 7° half-angle, the lepton is considered isolated if the sum of track and electromagnetic 
energy between the 7° half-angle cone and a 25° half-angle isolation cone does not exceed 
15 % of the lepton energy. 

For each individual centre-of-mass energy there is good agreement between the numbers of ex- 
pected background events and observed candidates after the preselection. Table 2 gives the 
number of preselected events summed over all centre-of-mass energies. Figure 1 shows the dis- 
tributions for background classes summed over all centre-of-mass energies and three arbitrarily 
scaled signal distributions (at a centre-of-mass energy of 206 GeV). The efficiencies of the pres- 
election vary on average between 39 % and 55 % for small decay widths and between 45 % and 
53 % for larger decay widths above Th = 100 GeV. 
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4.3 Likelihood analysis 

To consider the changing kinematic properties of the signal hypotheses in an optimal way, five 
different likelihood-based analyses for the signal and background discrimination were applied 
after the preselection. By a likelihood analysis we denote the combination of a set of likelihood 
input variables, a so-called likelihood, and the corresponding reference distributions of these 
variables. The reference distributions are filled with events of the specific classes for which the 
likelihood is calculated. The classes considered in this search are the two- and four-fermion 
backgrounds and the signal events. The two-photon events are negligible after the preselection. 
The search uses combinations of two likelihoods and three fixed signal mass ranges for unweighted 
reference histograms. 

To compare the kinematic properties of a selected data event to the hypothesis (Mh,Th) 
when evaluating the likelihood, one in principle has to fill weighted signal reference distributions 
for each hypothesis (Mh,Th). This will soon lead to an unmanageable technical effort, given the 
number of hypotheses scanned. Therefore a compromise was sought in which certain kinematic 
properties of the signal were emphasised and simultaneously the number of reference histograms 
kept small. This was achieved by filling unweighted signal reference histograms. For most of the 
(Mh,Th) hypotheses all signal masses were used for filling the reference histograms. This reflects 
the fact that for a very large decay width of the Higgs boson the possible values of kinematical 
variables are also smeared out over a large range. It was, however, found that the sensitivity of 
the likelihood selection (i.e. the median expected upper limits on a(Mn,Tu) x BR(H — > -Emis)) 
could be increased further for small widths below 50 GeV by filling reference histograms with 
signal masses from 50-80 GeV and from 80-120 GeV for intermediate and heavy Higgs boson 
masses respectively. A first likelihood was designed for a signal consisting of small masses (Mr 
< 80 GeV) or large masses and a very large width (Th > 110 GeV). In this likelihood input 
variables are used exploiting the characteristics of the dominant fraction of light masses in the 
signal mass distribution. However for signal masses above Mh = 80 GeV and small or moderate 
(i.e. below 110 GeV) decay widths, the contribution of large masses dominates the signal mass 
distribution. In this case the kinematics and topology of the signal events are determined by 
the higher masses close to the kinematic limit. A second likelihood is therefore built with input 
variables optimised for such signal characteristics. In the following the choice of the inputs for 
the two optimised likelihoods are presented. 

The first three input variables are used in both likelihoods (see Figure 3) . 

(1) (l+7>(M V is ^Mz))- 1 

V(Myis = Mz) is the probability of a kinematical x 2 fit of the jet four-vectors under the 
assumption that the invariant mass of the two jets is compatible with the Z boson mass. 
This variable depends only weakly on the Higgs mass. For events with non-converging fit 
the probability is set to zero. They therefore accumulate at a value of 1. 

(2) -logy 32 

The Durham algorithm groups two energy flow objects i and j into a jet as long as their 
separation in phase space yij = 2 x min(Ef,Ej ) x (1 — cos(#ij)/£y IS is smaller then the 
cut value y cu t. The number of jets in a event is predefined to be 2, 1/32 is the value of y cut 
where the two-jet topology of the event changes to a three-jet topology. Hence the negative 
logarithm of the so-called jet resolution parameter 1/32 is a measure for the jet topology 
being more two-jet like (large value of — log 1/32) or three-jet like (small value of — log 2/32). 

( 3 ) Pmis/v 7 * 

The transverse missing momentum p^is is one of the most prominent characteristics of 
signal-like events, but depends very much on the Higgs boson mass. For a heavy Higgs 
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boson produced close to the kinematic threshold almost at rest, the Z boson has almost no 
boost and decays into two more or less back-to-back jets. In this case the discriminating 
power of the variable is lost. 

The next three variables (see Figure 4) complete the first likelihood, which is used for all Higgs 
masses in the domain of very large width > 110 GeV or low Higgs masses < 80GeV. 

(4) ^ACOL 

The acolinearity angle <^>acol of the two-jet system is obtained by subtracting the three- 
dimensional angle between the reconstructed jet-axes from 180°. Events containing a low- 
mass Higgs boson exhibit on average a larger acolinearity than the background. 

(5) | cos(9*| 

The Gottfried- Jackson angle 9* , is defined as the angle between the flight direction of the 
Z boson in the laboratory frame and the direction of the decay products of the Z boson 
boosted into the Z boson rest-frame. The variable tends to have smaller values for the 
signal. 

(6) -logj/21 

The variable — log 2/21 is analogous to — log 1/32 and measures the compatibility of the event 
with a two-jet topology. Two-jet events tend to accumulate at small values of — log j/21 • 

The last three variables (see Figure 5) tune the second likelihood to become more sensitive for 
large Higgs boson masses and small to moderate widths. 

(7) Ef^/^S 

The variable Ej^f measures the energy of the most energetic of the two jets. This is on 
average higher for the four-fermion background, due to the boost of the W and Z pairs, 
whereas heavy Higgs bosons and a Z boson are produced at rest. 

(8) Hr ti 

This variable is the significance of the acoplanarity between the two jets, taking into account 
detector resolution and acceptance. The discrimination power is enhanced by weighting 
the acoplanarity with the average jet polar angle, since transverse jet directions are more 
precisely measured at large polar angles. Signal events tend to have a more significant 
acoplanarity and thus larger values of Rp ti than background. The precise definition of Rp ti 
can be found in the OPAL analyses of ZZ — ► qqi/77 events [25]. 

(9) (Mvis + M M is)/(M V is - M M is) 

This variable, described in [26], uses two strongly correlated quantities, the invariant miss- 
ing mass Mmis and the visible mass of the event Myis- Depending on the mass recon- 
struction accuracy it can have positive or negative values. The signal distribution of this 
variable is broader and accumulates at higher values than for the two- and four-fermion 
events, which are distributed more narrowly around the origin. 

From the two likelihoods and three ranges of signal masses filled in the reference histograms 
one has six analyses to search for the different hypotheses in Mh and Th- The study of the 
median expected <7(Mh,Fh) x BR(H — ► Emis) shows that five of these six are sufficient to have 
an optimally efficient analysis for each signal hypothesis characterised by Mh and Th (see Figure 
2) in the range studied. Likelihood 1 was not used with the reference distribution filled for the 
signal mass range of 80 to 120 GeV. Figure 6 a) to c) and g) to h) display examples for the 
likelihood distributions of all five analyses used. In the histograms the events selected at all 
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five centre-of-mass energies are added up, although each centre-of-mass energy was evaluated 
separately in the limit setting as explained in Section 5.1. The appropriate likelihood was cal- 
culated for each background, data and signal event. In case of a signal event it was added to 
the histogram with the weight defined in Equation 7. The number of expected signal events is 
normalised according to Equation 6. The use of different analyses gives rise to varying shapes of 
the likelihood distributions of the background. Also the various shapes of the signal likelihood 
distributions are visible for different Mh and Tr. Since the form of the likelihood distributions 
for signal and background can yield additional information in the limit calculation, only a loose 
cut is applied in the likelihood selection, requiring a signal likelihood larger than 0.2. 

4.4 Correction of background and signal efficiencies 

A correction is applied to the number of expected background events and the signal efficiencies 
due to noise in the detectors in the forward region which is not modelled by the Monte Carlo. 
The forward energy veto used in the preselection can accidentally be triggered by machine back- 
grounds. The correction factor is derived from the study of random beam crossings, and applied 
individually for each year of data taking. Random beam crossing events were recorded when no 
physics trigger was active. The fraction of events that fail the veto on activity in the forward 
region is below 3.4% for all runs analysed. The detailed breakdown of the fraction of accidentally 
vetoed events is given in the last column of Table 1. 

4.5 Systematic uncertainties 

A possible signal in the data would reveal itself by altering the shapes of the distributions of 
the discriminating variables. Thus a systematic deviation in the description of a reconstructed 
observable between Standard Model Monte Carlo and a data sample in which the signal is absent, 
could wrongly be attributed to the presence of a signal. 

The systematic uncertainties in the Monte Carlo description of the kinematic event variables 
are studied in two control samples at a centre-of-mass energy of 206 GeV. In the first control 
sample, called two-fermion control-sample in the following, radiative returns contributing to the 
qq(7) processes with photons detected at large angles are selected and the tagged ISR photon 
is removed from the event in Monte Carlo and data. This creates a qq-like topology with 
missing momentum at large angles. The second control sample, called four-fermion control- 
sample, is obtained by selecting W + W~ — ► qqlv events and removing the identified isolated 
lepton from the events in Monte Carlo and data. After this procedure control samples possess a 
topology very similar to signal events. For all kinematic variables x of the preselection and the 
likelihood selection the mean x and the width of the distribution (RMS) are compared between 
Monte Carlo and data, for the two-fermion and four-fermion control-samples. The observables 
in the two-fermion and four-fermion Monte Carlo are then modified separately according to 
x mc W = ( x mcP ~%c) x R rms^ +^data- Then all five likelihood selections are repeated 
separately and the relative change in the number of selected events compared to the unmodified 
case is taken as the systematic uncertainty. 

Since the analysis labelled Al in Table 3 is used over a large range of the search plane 
(see Figure 2), the systematic uncertainties on the background determined in this analysis Al 
are taken as an estimate for the background for all analyses. To determine the effect of the 
systematic uncertainty on the signal efficiencies, one has to take into account the fact that 
the kinematic properties of the signal depend on the assumed Higgs mass and decay width. 
Twelve representative hypotheses are studied with Mh chosen to be 20, 60 or 110 GeV and 
Th taking values of 5, 20, 70 and 200 GeV. For these hypotheses the signal Monte Carlo is 
modified according to the four-fermion correction factors, representing the dominant remaining 
background after the cut on the signal likelihood. The relative change in selected event weights 
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compared to the unmodified case is then taken as an estimate of the systematic uncertainty on 
the signal efficiency for a given hypothesis. The root-mean-square of all twelve hypotheses is 
applied as an (MH,rn)-independent estimate for the whole search area and for all centre-of-mass 
energies (see Table 4). 

The W pairs are very effectively reduced in the preselection by the isolated lepton veto. Due 
to the importance of this veto the uncertainty from the lepton isolation angle and the vetoed cone 
energy is studied in the following way. The half-cone angle of the outer cone is increased and 
decreased by two degrees, following the studies in [27], and the relative effect on the selection 
determined. Furthermore the cone energy is varied by 7.4 % and the analyses are repeated. 
The value of the cone energy rescaling is determined by the relative deviation of the mean of the 
measured energy of the lepton candidates in the inner cone between data and Monte Carlo in the 
W + W _ — > qqlv sample. For signal efficiencies an analogous study was performed at the twelve 
points described above. Both results for the relative change of the selection for the cone opening 
half-angle and cone energy variation are added in quadrature and the root-mean-square of the 12 
(Mh,Fh) hypotheses was taken to yield the total uncertainty associated with the isolated lepton 
veto (see Table 4). 

The theoretical prediction on the cross-section for the two- and four-fermion processes adds 
an uncertainty of 2% to the background uncertainty [28]. Finally, the uncertainty due to the 
limited Monte Carlo statistics is evaluated. 

Table 5 summarises the results of the studies. All uncertainties are assumed to be uncorre- 
cted and the individual contributions are added in quadrature to obtain the total systematic 
uncertainties on the background expectation and signal efficiency. The dominant systematic 
uncertainties on the signal efficiency arise from the description of the kinematic variables. The 
background expectation is more affected by the uncertainty in the isolated lepton veto, as the 
main contribution of the background stems from four-fermion processes. But the uncertainty 
associated with the description of the kinematic variables is of similar magnitude. The limits 
quoted in Section 5.1 were calculated including the uncertainties of Table 5. To estimate the 
extent to which the limits depend on the size of the systematic uncertainties, the limit calculation 
was repeated doubling the systematic uncertainties. A comparison of the limits with single and 
double systematic uncertainties, done at similar representative points as used for the systematic 
studies, showed that the excluded cross-sections typically decrease between a half and one and 
a half percent. A maximal reduction of 2.1 % was found. 

5 Results 

The results of the search using each of the five different likelihood selections, labelled A1-A5, after 
a cut on the likelihood larger than 0.2 are summarised in Table 6, which compares the numbers of 
observed candidates with the Standard Model background expectations. The data are compatible 
with the Standard Model background expectations. The remaining four-fermion background 
consists predominantly of W pairs, representing roughly three quarters of the background at 
energies above the Z pair threshold. Figure 7 shows examples for signal efficiencies. For small 
decay widths the dependence on the centre-of-mass energy is weak up to Mh ~ 80 GeV, and 
for large widths it is weak up to the kinematic limit. Because of the centre-of-mass energy 
dependence of the Bjorken cross-section of a Higgs boson with mass Mh, the lower centre-of- 
mass energies contribute more significantly to the sensitivity for lighter Higgs bosons (see e.g. 
in [8]). For very light Higgs bosons the efficiency is moderately reduced by the preselection cuts 
demanding a sizeable amount of missing energy. In the case of broader Higgs resonances with 
high mass, one observes a generally enhanced efficiency since the chance of selecting events from 
the low mass tail compensates the suppression due to the falling production cross-section of a 
heavy Higgs boson. 
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5.1 The upper limits on the production cross-section times branching ratio 

As no significant excess over the expected background is observed in the data, upper limits 
are calculated on the model-independent cross-section 0"(Mh,Fh) x BR(H — > -Emis) scaled to 
y/s = 206 GeV. As the likelihood distributions are only loosely cut, one can use not only 
the information from the integral number of selected events (Table 6) but also from the shape 
in a likelihood ratio [29] to set more sensitive upper limits. For each centre-of-mass energy 
separately, each bin with a signal likelihood larger than 0.2 in the distributions of expected 
signal, background and selected data is treated as a search channel. For each centre-of-mass 
energy the number of expected signal events is scaled to the total cross-section (Equation 6). As 
with the analysis described in [30], the likelihood distributions are given as discriminating input 
to a limit program [31]. A likelihood ratio is used to determine the signal confidence level, CLs, 
defined in [29, 31], which excludes the presence of a possible signal according to the modified 
frequentist approach [31]. Additionally the program calculates the median upper number of 
signal events that could be excluded at 95% confidence level (CL). This number is then scaled 
to the total cross-section at the centre-of-mass energy of 206 GeV for each (Mh,Th) hypothesis. 
The systematic uncertainties on the background expectations and signal selection efficiencies are 
included according a generalisation of the method described in [32]. 

A very fine scan of the (MH,rn)-plane was performed by simulating the spectra of Higgs 
bosons with a mass Mh from 1 to 120 GeV and widths Th starting at 1 GeV up to 3 TeV. The 
Higgs boson mass was simulated in steps of 5Mu = 1 GeV. Simulated values of Th are spaced 
in steps of 1 GeV up to 5 GeV. A spacing of 5Tn = 5 GeV is chosen from Th = 5 GeV to Th = 
750 GeV. Above this value steps of ST^ = 50 GeV are adopted up to the maximal Th of 3 TeV. 

Examples of the projections of the observed upper cross-section limits together with the 
median expected upper limits and the corresponding one and two standard deviation bands 
on the expected limits are displayed in Figure 8 for some choices of Th- Above a width of 
300 GeV the exclusion plots look quite similar to the example displayed in Figure 8i) because 
the excluded limits do not change very much. The observed limits for Th > 60 GeV are well 
contained in the one standard deviation bands on the expectations and generally do not exceed 
two standard deviations except in Figure 8 a) at Mh = 114 GeV. The discontinuities in the 
graphs correspond to changes in the analyses. As one can observe, below Th < 40 GeV the 
analysis are changed more often. Therefore the chance is higher that in a few bins there are 
statistical fluctuations in the selected data, that lead to a deviation of more than one standard 
deviation around the median. Also the data selected are highly correlated, as one can see for 
example in the upward fluctuation around Mh = 114 GeV visible in Figure 8 a)-c). All results 
for the observed upper limits on ct(Mh,Fh) x BR(H — > -Emis) are summarised in a contour plot 
(Figure 9) in the scanned (MH,rn)-plane. Above Th = 200 GeV the observed upper limits are 
in the range of 0.15 pb to 0.18 pb for all Mh and vary very little. For such large Th the recoil 
mass distribution of the Higgs tends to be more and more uniformly stretched out over the mass 
range explored. There is not much difference in the selection of signal events for a Higgs boson 
with e.g. a width of 400 GeV or 600 GeV in the considered range of Higgs masses. This prevents 
any specific discriminating kinematical properties from being assigned to the expected signal as 
signal masses of a broad kinematical range are selected with roughly equal probability. Therefore 
only one likelihood analysis is used in this part of the search area, selecting the same subset of 
data and background. Since the upper limit on the model-independent cross-section refers to a 
production cross-section at a centre-of-mass energy of 206 GeV, it must become independent of 
(Mh,Fh) for an extremely large Th. In this case the shape of the Higgs signal would just be 
a box, weighted with the production cross-section from 1 GeV to the kinematic limit of about 
115 GeV. The data are then compared to an approximately constant signal expectation. Hence 
the upper limit on the cross-section is approximately independent of the (Mh^h) hypothesis at 
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a value of roughly 0.16 pb. For resonances with a decay width smaller than 200 GeV there are 
regions where the limits are below 0.15 pb or even 0.1 pb for Mh between 60 and 74 GeV. In this 
mass range the number of data events selected is smaller than expected. Above Mh of 85 GeV 
the upper limits become larger than 0.2 pb and rise considerably for small widths below 40 GeV 
(see Figure 8 a-e). This is due to the fact that the Higgs mass approaches the kinematic limit 
and the likelihoods which rely on kinematical variables like p^is lose discrimination power. A 
maximal value of 0.57pb is observed for Mh of 114 GeV and Th of 1 GeV corresponding to a 
circa two-tr excess in the data. 

It should be kept in mind that no optimisation of the search has been performed for Th 
below 5 GeV. In the region of heavy Higgs boson mass > 105 GeV and small width a search 
using recoil mass spectra would be more sensitive. Therefore this region is more sensitively 
covered by searches that have been performed by the LEP experiments documented in [2]. 

5.2 Interpretation of the result in the stealthy Higgs scenario 

Interpreting the width Th of a Higgs boson according to Equation 5, and setting m p hi n to zero, 
it was possible to set limits on u in the stealthy Higgs scenario. A range from uj = 0.04 to 
to = 24.45 was probed. The excluded regions are shown in Figure 10 at 95 % confidence level 
(CL) in the w-Mh parameter space. To illustrate the Higgs boson width according to Equation 5, 
contours of fixed Th corresponding to a given mass Mh and coupling to are added to the plot. 
The maximum excluded invisible width is about Th = 400 GeV for Higgs boson masses < 35 GeV, 
decreasing slowly to Th = 115 GeV for Mh = 100 GeV. The minimal exclusion of uj = 0.04 is 
observed at Mh = 1 GeV and the maximal exclusion is u = 5.9 for Mh = 73 GeV. For uj between 
0.04 and 0.59 a Higgs mass from 1 to 103 GeV could be excluded. The maximal excluded Higgs 
mass was 103 GeV for width between 1 and 3 GeV, compared with the expected exclusion of 
106 GeV. 

The results presented in this study extend the previous decay-mode independent searches 
for new scalar bosons with the OPAL detector [5] to regions of larger couplings and higher 
Higgs boson masses. In [5] an interpretation within the stealthy Higgs model yielded a maximal 
excluded coupling uj for masses around 30 GeV, where uj was excluded up to uj = 2.7. That study 
excluded Higgs boson masses up to Mh = 81 GeV. It should be pointed out that the decay-mode 
independent searches also studied Higgs widths between 0.1 and 1 GeV and therefore cover the 
gap between searches within scenarios assuming a narrow decay width of the invisibly decaying 
Higgs boson [2] and the search presented in this paper up to Mh = 81 GeV. 

6 Conclusions 

A dedicated search was performed in the channel e + e~ — > HZ with Z — ► qq and the non- 
standard Model decay H — ► .Emis final state allowing for invisible decay widths of the Higgs 
boson from 1 GeV up to 3TeV. The data taken by the OPAL detector at LEP above the W 
pair threshold were analysed. No indication for a signal was found and upper limits were set on 
0-(M H ,r H ) x BR(H ->■ £mis)- The maximal upper limit is 0.57 pb at M H = 114GeV and T H = 
1 GeV. Over the scanned region of the (MH,rn)-plane upper limits are generally of the order of 
0.15 pb, especially for large values of Th > 400 GeV or Higgs boson masses < 85 GeV. 

The limits were interpreted in the stealthy Higgs scenario assuming the presence of a large 
number of massless singlet states. Limits were calculated on the coupling uj to a hidden scalar 
sector of the Higgs boson with a given mass Mh. A large part of the parameter plane kinemati- 
cally accessible with LEP 2 was excluded extending a previous exclusion published in [5] . Values 
for uj between 0.04 (Mh = 1 GeV) and 5.9 (Mh = 73 GeV) were excluded, and for certain values 
of uj Higgs boson masses are excluded up to Mh = 103 GeV. The possible non-detection of a 
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light Higgs boson at the LEP searches due to non-Standard Model invisible Higgs boson decays 
is therefore restricted to the case of extremely large decay widths > 400 GeV. 
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binned y/s 


nominal y/s (GeV) 


year 


int. luminosity (pb 1 ) 


accid. veto (%) 


> 180-186 


183 


1997 


40.0 


3.37 


> 186-193 


189 


1998 


199.8 


2.24 


> 193-198 


196 


1999 


70.4 


2.53 


> 198-203 


200 


1999 


112.0 


2.96 


> 203-209 


206 


2000 


206.9 


2.22 



Table 1: Breakdown of the analysed integrated data luminosities according to the centre-of-mass 
energies. The data was binned in five nominal centre-of-mass energies. The last column states 
the reduction of the signal efficiencies and expected background rates due to accidental triggering 
of the forward energy veto in the preselection, which is not modelled in the Monte Carlo. 



cut 


77 


qq( 7) 


4-fermion 


total SM 


data 


(l)-(5) 


48795 


15639 


4880 


69314 


74178 


(6) 


148 


10359 


1394 


11901 


11779 


(7) 


62 


9128 


1336 


10526 


10472 


(8) 


44 


4897 


1167 


6108 


6264 


(9) 


33 


1061 


964 


2058 


2116 


(10) 


18 


425 


895 


1338 


1387 


(11) 


18 


423 


879 


1320 


1368 


(12) 


4 


68 


820 


892 


899 


(13) 


4 


60 


441 


505 


498 



Table 2: Expected number of Standard Model background events after the preselection nor- 
malised to a data luminosity of 629.1 pb -1 . The total SM background after preselection is 
expected to be 505 ± 5(stat) ± 21(syst). The contributions of the different subclasses are bro- 
ken down in column two to four for the two-photon, two-fermion and four-fermion processes 
respectively. 



label 


likelihood 


analysis 

reference mass range (GeV) 


background uncertainty 
kinematic var. isol. lepton veto 


Al 


1 


1-120 


2.4% 


2.4% 


A2 


2 


1-120 


1.6% 


2.3% 


A3 


1 


50-80 


1.0% 


2.5% 


A4 


2 


50-80 


1.6% 


2.6% 


A5 


2 


80-120 


1.1% 


1.5% 


choice for uncertainty 


2.4% 


2.4% 



Table 3: Results of the study of systematic uncertainties of the expected background for the 
five kinds of analyses, labelled A1-A5, used in the search (see Figure 2) at a centre-of-mass 
energy of 206 GeV. Since the analysis labelled Al covers the largest part of the search area, its 
uncertainty was chosen as representative uncertainty on the background due to the uncertainty 
in the kinematic variables and the isolated lepton veto at all centre-of-mass energies. 
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signal hypothesis 
M H (GeV) T H (GeV) 


efficiency uncertainty 
kinematic var. isol. lepton veto 


zO 


5 


0.6 7o 


0.6 7o 


on 


on 


U.I A) 


U. 1 A) 


20 


70 


n 3 % 

u.o /u 


7% 

U. ( /u 


20 


200 


0.1% 


0.7% 


60 


5 


0.7% 


0.8% 


60 


20 


0.7% 


0.8% 


60 


70 


0.2% 


0.8% 


60 


200 


0.3% 


0.7% 


110 


5 


5.5% 


0.7% 


110 


20 


2.9% 


0.8% 


110 


70 


1.3% 


0.8% 


110 


200 


0.1% 


0.8% 


all M u 


and Th 


1.9% 


0.7% 



Table 4: Results of the study of systematic uncertainties in twelve representative (Mu,Tu)-points 
at a centre-of-mass energy of 206 GeV. For each source the root-mean-square of the individual 
uncertainties in the twelve points was taken to get an (Mn,Tn) independent estimate of the 
uncertainty at all centre-of-mass energies. 



source 


background uncertainty 


efficiency uncertainty 


kinematic variables 


2.4% 


1.9% 


isolated lepton veto 


2.4% 


0.7% 


limited MC statistics 


1.0% 


0.2% 


prediction 2- and 4-f cross-sect. 


2.0% 




total uncertainty 


4.1% 


2.0% 



Table 5: Results of the study of systematic uncertainties of the background for the five analyses 
(see Table 3) and of the signal efficiencies in twelve representative (Mn,Tu)-points (see Table 4) 
at a centre-of-mass energy of 206 GeV. The total uncertainty on background expectation and 
signal efficiency is applied at all centre-of-mass energies and for all (Mn,Tu) hypotheses. 



label 


reference masses 


likelihood 


2-fermion 


4-fermion 


total SM 


data 


Al 


1-120 GeV 


type 1 


11 


374 


385 ± 4 ± 16 


369 


A2 


1-120 GeV 


type 2 


3 


378 


381 ± 4 ± 16 


370 


A3 


50-80 GeV 


type 1 


5 


315 


320 ± 3 ± 13 


305 


A4 


50-80 GeV 


type 2 


2 


315 


317 ± 3 ± 13 


310 


A5 


80-120 GeV 


type 2 


8 


247 


255 ± 3 ± 11 


253 



Table 6: The likelihood selection of events with a signal likelihood exceeding 0.2 according 
to the different search strategies. The individual contributions to the total Standard Model 
background of two-fermion and four-fermion background is broken down in the second and third 
column respectively. For the total Standard Model background the statistical and the systematic 
uncertainty is also given. 
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Figure 1: Distribution of the preselection variables after the preselection cuts (l)-(5). All classes 
of Standard Model background and data are added for all analysed centre-of-mass energies. The 
distributions of three arbitrarily scaled signal hypotheses at = 206 GeV are displayed as open 
histograms. 
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used likelihood No. 1 for this hypothesis 

used likelihood No. 2 for this hypothesis 
used reference mass range 1-120 GeV for this hypothesis 
used reference mass range 50-80 GeV for this hypothesis 
used reference mass range 80-120 GeV for this hypothesis 



Figure 2: In total five analyses were used to cover the plane of hypothetical Higgs mass and decay 
width pairings. The analyses differ in whether the first or second likelihood was used (denoted 
by the number in the cell) and what signal masses where used in filling the reference histograms 
(depicted by the shading of the cell). The pattern resulted from an optimisation starting with 
= 5 GeV up to 50 GeV. Below Th = 5 GeV the pattern was simply continued and not 
optimised anymore. Above 50 GeV a simple continuation of the pattern was found and proved 
to be sufficiently sensitive. 
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Variables used in likelihoods 1 and 2 

• OPAL 

I I 2-fermion 
4-fermion 



(M H ,r H )=(110,10) (GeV) 
-■ (M H ,r H )=( 80,15) (GeV) 
— (M H ,r H )=( 40,20) (GeV) 




(1+P(M VIS =M 7 )) 




Figure 3: Distributions of the likelihood variables. All classes of Standard Model background 
and data are added for all centre-of-mass energies analysed. The distributions of three arbitrary 
scaled signal examples at yfs = 206 GeV are displayed as open histograms. The variables shown 
contribute to likelihood 1 and 2 as they exploit general properties of the signal signature. 
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lcos9*l " lo §ioyi2 

Figure 4: Distributions of the likelihood variables. All classes of Standard Model background 
and data are added for all centre-of-mass energies analysed. The distributions of three arbitrary 
scaled signal examples at y/s = 206 GeV are displayed as open histograms. The variables shown 
are combined with the ones of Figure 3 to construct the likelihood 1 used in a general search at 
different Mh and Th- 
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Variables used in likelihood 2 

• OPAL 
I I 2-fermion 
I I 4-fermion 

(M H ,r H )=(110,10) (GeV) 
-- (M H ,r H )=( 80,15) (GeV) 
— (M H ,r H )=( 40,20) (GeV) 





j III! rrrr+4— WW II II 

-50 -40 -30 -20 -10 10 20 30 40 50 



(M VIS+ M MIS )/(M VIS -M MIS ) 



Figure 5: Distributions of the likelihood variables. All classes of Standard Model background 
and data are added for all centre-of-mass energies analysed. The distributions of three arbitrary 
scaled signal examples at yfs = 206 GeV are displayed as open histograms. The variables shown 
have a larger discrimination power for a heavier Higgs boson and contribute with the variables 
of Figure 3 to the second likelihood. 
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Figure 6: Examples of some of the likelihood selections. Figure 6 e), f), h) corresponds to analysis 
Al (as labelled in Table 3 and 6), a) to A3, b) to A4 and c), d) to A5. The OPAL data and the 
expected 2-fermion and 4-fermion background are added for all analysed centre-of-mass energies. 
The signal hypothesis in the hatched histograms is normalised to the number of expected signal 
events and added to the background. The open histograms display the shapes of scaled signal 
distributions. 
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Figure 7: Examples for the selection efficiency after a cut on the signal likelihood greater than 
0.2 versus the Higgs mass Mh as function of the assumed decay width Th at the different y/s. 
The error is the binomial error on the selected event weights and smaller than the markers. Lines 
are added to guide the eye.For a large Th signal hypothesis, the kinematic distributions of events 
and the distribution of weights assigned to this events are broader. Therefore it is more likely 
to select a larger fraction of the event weights leading to a higher and more uniform efficiency. 
For a smaller widths Th the efficiency to detect relatively heavy and a more Standard Model like 
Higgs boson is more restricted by the available \/s. 
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Figure 8: The model independent upper limits at 95 % CL on the production cross-section times 
branching ratio, <7 pro d X BR mv , scaled to a centre-of-mass energy of 206 GeV for Higgs mass Mr 
and some examples of the Higgs decay width Th. The discontinuities in the limits reflect the 
changes in the analysis used at this mass (see Figure 2). 
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Figure 9: The exclusion contours at 95 % CL on the model independent production cross-section 
times branching ratio, a pro d x BR mv , for Higgs boson mass Mh and the Higgs boson decay width 
Th up to 3 TeV (note a change in logarithmic scale below Th = 5 GeV for better visibility). Solid 
and dashed lines delimit areas of excluded upper limits. Cross-sections times branching ratio 
between 0.07 pb and 0.57 pb are excluded with the OPAL data above yfs = 183 GeV. 
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Figure 10: The exclusion contours at 95% CL on the Higgs mass Mh and the Higgs-phion 
coupling u of the stealthy Higgs model. The values of u are related to the decay width Th via 
Th(Mh) = Fsm(-^h) + 32*7]- Mh ' * n ^ e case °^ mass ^ ess phions (see Equation 5). Contours of 
fixed Th are also shown in the plot as dashed lines. 



30 



EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



OPAL PR417 
CERN-EP-2006-026 
August 3, 2006 



Search for Invisibly Decaying Higgs Bosons 

with Large Decay Width Using the 
OPAL Detector at LEP 

The OPAL Collaboration 



Abstract 

This paper describes a topological search for an invisibly decaying Higgs boson, H, produced via 
the Bjorken process (e + e~ — ► HZ). The analysis is based on data recorded using the OPAL de- 
tector at LEP at centre-of-mass energies from 183 to 209 GeV corresponding to a total integrated 
luminosity of 629 pb -1 . In the analysis only hadronic decays of the Z boson are considered. A 
scan over Higgs boson masses from 1 to 120 GeV and decay widths from 1 to 3000 GeV revealed 
no indication for a signal in the data. From a likelihood ratio of expected signal and Standard 
Model background we determine upper limits on cross-section times branching ratio to an invis- 
ible final state. For moderate Higgs boson decay widths, these range from about 0.07 pb (Mg = 
60 GeV) to 0.57 pb (Mh = 114 GeV). For decay widths above 200 GeV the upper limits are of the 
order of 0.15 pb. The results can be interpreted in general scenarios predicting a large invisible 
decay width of the Higgs boson. As an example we interpret the results in the so-called stealthy 
Higgs scenario. The limits from this analysis exclude a large part of the parameter range of this 
scenario experimentally accessible at LEP 2. 



(Submitted to Eur. Phys. J.) 



The OPAL Collaboration 



G. Abbiendi 2 , C.Ainsley 5 , P.F. Akesson 7 , G.Alexander 21 , G. Anagnostou 1 , K.J.Anderson 8 , 
S.Asai 22 , D.Axen 26 , I.Bailey 25 , E. Barberio 7 ' p , T. Barillari 31 , R.J. Barlow 15 , R.J.Batley 5 , 

P.Bechtle 24 , T.Behnke 24 , K.W.Bell 19 , P.J.Bell 1 , G.Bella 21 , A. Bellerive 6 , G.Benelli 4 , 
S.Bethke 31 , O.Biebel 30 , O.Boeriu 9 , P. Bock 10 , M. Boutemeur 30 , S.Braibant 2 , R.M.Brown 19 , 

H. J. Burckhart 7 , S.Campana 4 , P. Capiluppi 2 , R.K.Carnegie 6 , A.A.Carter 12 , J.R.Carter 5 , 
C.Y.Chang 16 , D.G.Charlton 1 , C.Ciocca 2 , A. Csilling 28 , M.Cuffiani 2 , S.Dado 20 , A.De Roeck 7 , 

E.A.De Wolf 7 ' s , K.Desch 24 , B.Dienes 29 , J.Dubbert 30 , E. Duchovni 23 , G.Duckeck 30 , 
I.P.Duerdoth 15 , E.Etzion 21 , F.Fabbri 2 , P.Ferrari 7 , F.Fiedler 30 , I. Fleck 9 , M.Ford 15 , A.Frey 7 , 
P.Gagnon 11 , J.W.Gary 4 , C. Geich-Gimbel 3 , G. Giacomelli 2 , P. Giacomelli 2 , M.Giunta 4 , 
J.Goldberg 20 , E. Gross 23 , J. Grunhaus 21 , M.Gruwe 7 , A.Gupta 8 , C.Hajdu 28 , M.Hamann 24 , 
G.G.Hanson 4 , A.Harel 20 , M. Hauschild 7 , C.M.Hawkes 1 , R.Hawkings 7 , G.Herten 9 , 
R.D.Heuer 24 , J.C.Hill 5 , D. Horvath 28 ' c , P. Igo-Kemenes 10 , K.Ishii 22 , H. Jeremie 17 , 
P. Jovanovic 1 , T.R. Junk 6 '*, J. Kanzaki 22 '", D.Karlen 25 , K.Kawagoe 22 , T.Kawamoto 22 , 
R.K.Keeler 25 , R.G. Kellogg 16 , B.W.Kennedy 19 , S.Kluth 31 , T. Kobayashi 22 , M.Kobel 3 '*, 
S.Komamiya 22 , T.Kramer 24 , A. Krasznahorkay Jr. 29 - 6 , P. Krieger 6 -* , J. von Krogh 10 , T.Kuhl 24 , 
M.Kupper 23 , G.D. Lafferty 15 , H. Landsman 20 , D.Lanske 13 , D. Lellouch 23 , J.Letts , 
L.Levinson 23 , J.Lillich 9 , S.L.Lloyd 12 , F.K. Loebinger 15 , J.Lu 26 - b , A.Ludwig 3 '*, J.Ludwig 9 , 
W.Mader 3 '*, S. Marcellini 2 , A.J.Martin 12 , T.Mashimo 22 , P. Mattig™, J.McKenna 26 , 
R.A.McPherson 25 , F.Meijers 7 , W.Menges 24 , F.S. Merritt 8 , H.Mes 6 ' a , N.Meyer 24 , 
A.Michelini 2 , S.Mihara 22 , G. Mikenberg 23 , D.J.Miller 14 , W.Mohr 9 , T.Mori 22 , A. Mutter 9 , 
K.Nagai 12 , I. Nakamura 22 -", H.Nanjo 22 , H.A.Neal 32 , R.Nisius 31 , S W. O'Neale 1 -*, A. Oh 7 , 
M.J.Oreglia 8 , S.Orito 22 '*, C.Pahl 31 , G. Pasztor 4,9 , J.R. Pater 15 , J.E. Pilcher 8 , J. Pinfold 27 , 
D.E. Plane 7 , O.Pooth 13 , M. Przybycieh 7 ' n , A.Quadt 31 , K. Rabbertz 7 ' r , C.Rembser 7 , 
PRenkel 23 , J.M.Roney 25 , A.M.Rossi 2 , Y.Rozen 20 , K.Runge 9 , K.Sachs 6 , T.Saeki 22 , 

E.K.G.Sarkisyan 7 -^, A.D. Schaile 30 , O.Schaile 30 , P. Scharff-Hansen 7 , J.Schieck 31 , 
T. Schorner-Sadenius 7 ^, M.Schroder 7 , M.Schumacher 3 , R. Seuster 13 '^, T.G. Shears 7 '^ 1 , 
B.C.Shen 4 , P.Sherwood 14 , A.Skuja 16 , A.M.Smith 7 , R.Sobie 25 , S. Soldner-Rembold 15 , 
F.Spano 8 ^, A.Stahl 13 , D.Strom 18 , R. Strohmer 30 , S.Tarem 20 , M. Tasevsky 7 ' d , R.Teuscher 8 , 
M.A.Thomson 5 , E. Torrence 18 , D.Toya 22 , P.Tran 4 , I. Trigger 7 '"', Z. Trocsanyi 29 ' 6 , E.Tsur 21 , 
M.F.Turner- Watson 1 , I.Ueda 22 , B. Ujvari 29 ' 6 , C.F. Vollmer 30 , P.Vannerem 9 , R. Vertesi 29,e , 
M.Verzocchi 16 , H.Voss 7 '"?, J. Vossebeld 7 -' 1 , C.P.Ward 5 , D.R.Ward 5 , P.M. Watkins 1 , 
A.T.Watson 1 , N.K.Watson 1 , P.S.Wells 7 , T.Wengler 7 , N.Wermes 3 , G.W. Wilson 15 ' fc , 
J.A.Wilson 1 , G.Wolf 23 , T.R.Wyatt 15 , S. Yamashita 22 , D.Zer-Zion 4 , L. Zivkovic 20 



School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 
2 Dipartimento di Fisica dell' Universita di Bologna and INFN, 1-40126 Bologna, Italy 
3 Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 
4 Department of Physics, University of California, Riverside CA 92521, USA 
5 Cavendish Laboratory, Cambridge CB3 OHE, UK 

6 Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ottawa, 
Ontario K1S 5B6, Canada 

7 CERN, European Organisation for Nuclear Research, CH-1211 Geneva 23, Switzerland 
8 Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago IL 60637, 
USA 

9 Fakultat fur Physik, Albert-Ludwigs-Universitat Freiburg, D-79104 Freiburg, Germany 
10 Physikalisches Institut, Universitat Heidelberg, D-69120 Heidelberg, Germany 
n Indiana University, Department of Physics, Bloomington IN 47405, USA 



1 



12 Queen Mary and Westfleld College, University of London, London El 4NS, UK 
13 Technische Hochschule Aachen, III Physikalisches Institut, Sommerfeldstrasse 26-28, D-52056 
Aachen, Germany 

14 University College London, London WC1E 6BT, UK 

15 Department of Physics, Schuster Laboratory, The University, Manchester M13 9PL, UK 
16 Department of Physics, University of Maryland, College Park, MD 20742, USA 
17 Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec H3C 3J7, Canada 
18 University of Oregon, Department of Physics, Eugene OR 97403, USA 
19 CCLRC Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 
20 Department of Physics, Technion- Israel Institute of Technology, Haifa 32000, Israel 
21 Department of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
22 International Centre for Elementary Particle Physics and Department of Physics, University 
of Tokyo, Tokyo 113-0033, and Kobe University, Kobe 657-8501, Japan 
23 Particle Physics Department, Weizmann Institute of Science, Rehovot 76100, Israel 
24 Universitat Hamburg/DESY, Institut fur Experimentalphysik, Notkestrasse 85, D-22607 Ham- 
burg, Germany 

25 University of Victoria, Department of Physics, P O Box 3055, Victoria BC V8W 3P6, Canada 
26 University of British Columbia, Department of Physics, Vancouver BC V6T 1Z1, Canada 
27 University of Alberta, Department of Physics, Edmonton AB T6G 2J1, Canada 
28 Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P O Box 49, Hungary 
29 Institute of Nuclear Research, H-4001 Debrecen, P O Box 51, Hungary 

30 Ludwig-Maximilians-Universitat Miinchen, Sektion Physik, Am Coulombwall 1, D-85748 Garch- 
ing, Germany 

31 Max-Planck-Institute fur Physik, Fohringer Ring 6, D-80805 Miinchen, Germany 
32 Yale University, Department of Physics, New Haven, CT 06520, USA 



a and at TRIUMF, Vancouver, Canada V6T 2A3 
b now at University of Alberta 

c and Institute of Nuclear Research, Debrecen, Hungary 

d now at Institute of Physics, Academy of Sciences of the Czech Republic 18221 Prague, Czech 
Republic 

6 and Department of Experimental Physics, University of Debrecen, Hungary 
f and MPI Miinchen 

9 and Research Institute for Particle and Nuclear Physics, Budapest, Hungary 
h now at University of Liverpool, Dept of Physics, Liverpool L69 3BX, U.K. 

* now at Dept. Physics, University of Illinois at Urbana-Champaign, U.S.A. 
J and Manchester University Manchester, M13 9PL, United Kingdom 

k now at University of Kansas, Dept of Physics and Astronomy, Lawrence, KS 66045, U.S.A. 

1 now at University of Toronto, Dept of Physics, Toronto, Canada 

m current address Bergische Universitat, Wuppertal, Germany 

n now at University of Mining and Metallurgy, Cracow, Poland 

° now at University of California, San Diego, U.S.A. 

p now at The University of Melbourne, Victoria, Australia 

q now at IPHE Universite de Lausanne, CH-1015 Lausanne, Switzerland 

r now at IEKP Universitat Karlsruhe, Germany 

s now at University of Antwerpen, Physics Department, B-2610 Antwerpen, Belgium; supported 
by Interuniversity Attraction Poles Programme - Belgian Science Policy 

* now at Technische Universitat, Dresden, Germany 

u and High Energy Accelerator Research Organisation (KEK), Tsukuba, Ibaraki, Japan 



2 



v now at University of Pennsylvania, Philadelphia, Pennsylvania, USA 
w now at TRIUMF, Vancouver, Canada 
x now at DESY Zeuthen 

* now at CERN 
2 now at DESY 

* Deceased 



3 



1 Introduction 



An intense search for the Higgs boson was undertaken by all of the four LEP experiments in 
various Standard Model and non-Standard Model search channels. Searches for the Standard 
Model Higgs boson, exploiting the prediction for its decay modes and also searches for invisible 
Higgs boson decays as predicted by various extensions of the Standard Model, excluded Higgs 
masses up to 114.4 GeV [1, 2]. These latter searches assumed a rather small invisible decay width 
comparable to the predicted Standard Model decay width for a light Higgs boson and well below 
the experimentally achievable mass resolution of about 3 to 5 GeV. 

Recent theories that postulate the existence of additional spatial dimensions offer a new 
possibility for invisible Higgs decays [3]. In such theories the Planck mass is lowered to the TeV 
range and a rich spectrum of new particles appears, like graviscalars in the case of flat extra 
dimensions. Hence the Higgs boson can mix with the graviscalars, which leads to a missing 
energy signature in the detector [3]. This mixing can result in a large invisible decay width of 
the Higgs boson, depending on the model parameters, and would therefore alter the Standard 
Model branching ratios. As a consequence of the broadening of the Higgs resonance in the 
recoil mass spectrum, the signal-to-background ratio can deteriorate significantly. In a worst 
case scenario, searches optimised under the assumption of a narrow Higgs resonance might have 
missed the detection of a kinematically accessible Higgs boson at LEP. 

This paper describes a search for the Higgs boson, H, which imposes no constraints on the 
total decay width. The Higgs boson is assumed to be produced in association with a Z boson 
via the Bjorken process, e + e~ — > HZ, where the Z is required to decay hadronically and the 
invisible Higgs boson is detected as missing energy Emis in the event. The results are presented 
in a model-independent way in terms of limits on the Bjorken production cross-section times 
branching ratio, a(Mn,Tn) x BR(H — > Emis), at a centre-of-mass energy of 206 GeV, where 
Th is the Breit-Wigner width of the Higgs boson. A simple model extending the Standard 
Model with additional SU(3)cxSU(2)lxU(1)y singlet fields which interact strongly with the 
Higgs boson ("stealthy Higgs scenario" [4]) is chosen as an example for the interpretation of the 
result. This interaction gives rise to a large invisible decay width of the Higgs boson. This 
dedicated search expands on the previous decay-mode-independent search [5] carried out by the 
OPAL Collaboration which reported for the first time limits on the HZ production cross-section, 
interpreted in the stealthy Higgs model. 

The paper is organised as follows. Section 2 introduces the stealthy Higgs scenario. Section 
3 gives details about the modelling of signal and background. Section 4 describes the search and 
the results are interpreted in Section 5. We summarise the results in Section 6. 

2 The stealthy Higgs scenario 

In general renormalisable theories there might be other fundamental scalars, in addition to the 
Standard Model Higgs boson, that do not interact with normal matter. To investigate the influ- 
ence of a hidden scalar sector on the Higgs observables the stealthy Higgs scenario conjectures 
the existence of additional SU(3)c><SU(2)lxU(1)y singlet fields called phions. Radiative cor- 
rections to weak processes are not sensitive to the presence of singlets in the theory because no 
Feynman graphs containing singlets appear at the one-loop level. Since effects at the two-loop 
level are below the experimental precision, the presence of a singlet sector is not ruled out by any 
LEP 1 precision data [4]. These phions would not interact via the strong or electro- weak forces, 
but couple only to the Higgs boson [4], thus offering invisible decay modes to the Higgs. The 
width of the Higgs resonance can become large if either the number of such singlets, N, or the 
coupling, lo, is large, thus leading to a broad mass spectrum recoiling against the reconstructed 
Z boson. 
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The Lagrangian of the scalar sector in this model contains only four additional parameters 
compared to the Standard Model. Not listing the unchanged couplings of the Higgs boson to the 
fermions, the scalar part of the Lagrangian is given by 

-^scalar — ^Higgs ~i~ ^phion -^interaction (1) 

£ffi g gs = -d^d^-x^cp- y) 2 (2) 

£ P hion = ~ 9^0 - ^m 2 phion f - (iff (3) 

^ -2 t 

^interaction = ~ 7= 4> <P- (4) 

2v N 



The term £ffiggs describes the usual Standard Model Higgs doublet 4> acquiring the Standard 
Model vacuum expectation value, v, and having its self-coupling A. In the free Lagrangian of 
scalar singlets, -Cphion, the singlets with mass m p hi n are denoted as the 0(./V)-symmetric multiplet 
(p. The phions also have a self-coupling k, which is fixed at k(2Mz) = 0, to allow for the widest 
parameter range of the model. The self coupling term entering loop calculations is suppressed like 
1/N. The interaction term between the Higgs and the additional phions, interaction, leads to the 
phenomenological consequence of invisible Higgs decays because the Higgs boson couples to the 
phions independently of their mass. The strength of the coupling is instead proportional to the 
coupling constant lo, which is a free parameter of the model. Even though the vacuum-induced 
mass term of the phions after the symmetry-breaking is suppressed like 1/y/N [4], the phions 
occur in loop corrections to the Higgs boson propagator and therefore affect the resonance width 
of the Higgs boson. An analytic expression [6] for the change in the Higgs width compared to 
the Standard Model decay width, Tsm, can be found in the limit N — ► oo, when neglecting the 
self-coupling of the phions as a small effect: 

2 2 

r H (M H ) = r SM (M H ) + x y/l - Am^JMl (5) 

The cross-section for the Bjorken process can be calculated from Equations 9 and 10 of refer- 
ence [4]. Using the parametrisation for the invisible decay width (Equation 1 and 2 in [7]) one 
can express the total cross-section for the production and invisible decay by 

dSl ^e-^ZH)^) 7l{{M 2 { _ Sl / +SlT 2 i y (6) 

Here si denotes the invariant mass squared of the invisible decay products of the Higgs boson. 
The production rate of these invisible masses is given by the Standard Model cross-section 1 
cr( e + e -^zH)( s i") f° r a Higgs boson of mass y/sj. Hence the Standard Model cross-section com- 
pletely determines the dependence of the total cross-section on the centre-of-mass energy (see e.g. 
in [8]). Therefore the total cross-section goes rapidly to zero for Higgs boson masses above the 
kinematic limit. The effect of the convolution with the Breit-Wigner-like function is a broadening 
of the resonance in the recoil mass spectrum and hence a dilution of the signal-to-background 
ratio. In extreme cases of large invisible decay width one could expect the Higgs recoil mass spec- 
trum to mimic the background. In such extreme cases even a light and kinematically accessible 
Higgs boson might have escaped detection at LEP. 

In Section 5.2 we derive limits on the stealthy Higgs model under the assumption of m pn ion = 
0. By simulating signal spectra for different Higgs boson masses Mh and widths Th, we set limits 
in the cj-Mh plane in the large N limit. 

x By choosing u) > one can prevent the phions from acquiring a non-zero vacuum expectation value and avoid 
a Higgs-phion mixing due to a non-diagonal mass matrix. In case of non-zero mixing, the couplings of the lightest 
scalar to the gauge boson would decrease proportional to the cosine of the mixing angle. As a consequence the 
cross-section of the Bjorken process would be lowered. 
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3 Data sets and Monte Carlo samples 



3.1 The OPAL detector and event reconstruction 

The OPAL detector [9], operated between 1989 and 2000 at LEP, had nearly complete solid 
angle 3 coverage and excellent hermeticity. The innermost detector of the central tracking was a 
high-resolution silicon microstrip vertex detector [10] which lay immediately outside the beam 
pipe. The silicon microvertex detector was surrounded by a high precision vertex drift chamber, 
a large volume jet chamber, and z-chambers which measured the z coordinates of tracks, all in a 
uniform 0.435 T axial magnetic field. A lead-glass electromagnetic calorimeter with presampler 
was located outside the magnet coil. In combination with the forward calorimeters, a forward 
ring of lead-scintillator modules (the "gamma catcher"), a forward scintillating tile counter [9, 11], 
and the silicon-tungsten luminometer [12], the calorimeters provided a geometrical acceptance 
down to 25 mrad from the beam direction. The silicon-tungsten luminometer served to measure 
the integrated luminosity using small angle Bhabha scattering events [13]. The magnet return 
yoke was instrumented with streamer tubes and thin gap chambers for hadron calorimetry and 
is surrounded by several layers of muon chambers. 

The analysis is based on data collected with the OPAL detector at LEP 2 from 1997 to 
2000 at centre-of-mass energies between 183 and 209 GeV. The integrated luminosity analysed 
is 629.1 pb" 1 . To compare with the Standard Model Monte Carlo the data are binned in five 
nominal centre-of-mass-energy points, corresponding to the energies at which the Monte Carlo 
is produced, as detailed in Table 1. 

A fast online filtering algorithm classifies the events as multi-hadronic. Events are re- 
constructed from tracks and energy deposits ("clusters") in the electromagnetic and hadronic 
calorimeters. All tracks and energy clusters satisfying quality requests similar to those de- 
scribed in [14] are associated to form "energy flow objects". The measured energies are corrected 
for double counting of energy in the tracking chambers and calorimeters by the algorithm de- 
scribed in [14]. Global event variables, such as transverse momentum and visible mass, are 
then reconstructed from these objects and all events are forced into a two-jet topology using the 
Durham algorithm [15]. 

3.2 Signal and background modelling 

To determine the detection efficiency for a signal from an invisibly decaying Higgs boson and the 
amount of expected background from Standard Model processes, several Monte Carlo samples 
are used. Signal events for a hypothetical Higgs boson mass Mh decaying with arbitrary broad 
width Th are simulated by reweighting invisibly decaying events of type H — ► XiXi- The mass 
of neutralinos x? 1S chosen such that the Higgs boson with mass vrii can decay into a pair of 
neutralinos, which leave the detector without being detected. These Higgs bosons with decays 
into 'invisible' particles are generated with masses from 1 GeV to 120 GeV with the HZHA [16] 
generator. The HZHA events are generated assuming the Standard Model production cross- 
section <7( e + e -^ ZH ) for the Higgs boson. The test masses nii are spaced in steps of 1 GeV. The 
spacing of the test masses is chosen such that they are not resolved by the detector in the signal 
yielded after a reweighting procedure described in the following. From Equation 6 one extracts 
the event weights Wi(mi; Mh,Th) for a mass point nii contributing to the search for a Higgs 
boson of mass Mh and total decay width Th- The total decay width Th is defined as the sum of 

3 OPAL used a right-handed coordinate system. The z axis pointed along the direction of the electron beam 
and the x axis was horizontal pointing towards the centre of the LEP ring. The polar angle 9 was measured with 
respect to the z axis, the azimuthal angle <j> with respect to the x axis. 
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Standard Model width and invisible width . 



Wi(mf, Mh, Th) = 120Ge v l^T—, ( ? ) 

2^m J= lGeV dnij \ m 3> 



drrii 



da ( ff(e+e-^ZH)(m t ) 2mfrg v 

dm^ 7r((M 2 -m 2 ) 2 + m 2 r 2 )" l8j 

The Standard Model cross-section 0"( C + C -^ ZH ) for the Bjorken production process in Equation 8 
propagates the centre-of-mass energy dependence of the total cross-section into the weights. The 
unweighted signal Monte Carlo samples contain 2000 events per mass point m^. In the reweighted 
signal Monte Carlo sample all test masses contribute according to their weight. The reweighted 
masses Mh range from 1 to 120 GeV spaced in steps of 1 GeV. The smallest width simulated 
by this procedure is a Th of 1 GeV and the largest a Th of 3TeV. The detection efficiency for 
a Higgs boson with Mh and Th is estimated by the sum of selected event weights assuming 
binomial errors. 

The classes of Standard Model background processes considered are two-photon 2 , two- and 
four-fermion processes. For simulation of background processes the following generators are 
used: KK2F [17] and PYTHIA [18] (qq(7)), GRC4F [19] (four-fermion processes), PHOJET [20], 
HERWIG [21], Vermaseren [22] (hadronic and leptonic two-photon processes). For Monte Carlo 
generators other than HERWIG, the hadronisation is done using JETSET 7.4 [18]. The integrated 
luminosity of the main background Monte Carlo samples is at least 15 times the statistics of the 
data for the two-fermion background, 24 times for the four-fermion background and 30 times for 
the two-photon background. The Monte Carlo events are passed through a detailed simulation 
of the OPAL detector [23] and are reconstructed using the same algorithms as for the real data. 



4 Search for e + e — > HZ with Z — > qq and H — > .Emis final state 

The event selection is intended to be efficient for the complete range of possible Higgs masses 
Mh and corresponding decay widths Th studied in this search. The preselection cuts remain 
relatively loose and are intended to accumulate signal like event topologies in the data. The 
final discrimination between signal and background is done by a likelihood-based selection. The 
optimised likelihood selection has to account for the fact that the kinematical properties of the 
signal change considerably over the range of masses and width hypotheses considered. 



4.1 Event topologies 

The signal signature is generally characterised by an acoplanar two-jet system from the Z boson 
decay. We use the term 'acoplanar' for jet pairs if the two jet axes and the beam axis are 
not consistent with lying in a single plane. The decay products of the Z boson are preferentially 
emitted into the central part of the detector, recoiling against the invisibly decaying Higgs boson. 
This is because, in contrast to the irreducible background of ZZ — > qq vv which is produced 
with an angular dependence of the differential cross-section proportional to cos 2 9, the Bjorken 
process is proportional to sin 2 6. The Higgs boson decay leads to a large missing momentum and 
a significant amount of missing energy. In two-photon processes, where the incoming electron 
and positron are scattered at low angles, usually one or both of the electrons remain undetected. 
Events of this type have large missing momentum with the missing momentum vector, puis, 

2 Two-photon interactions occur when an electron and a positron at high energies and in close proximity emit 
a pair of photons which interact via the electromagnetic force to generate a fermion pair. 
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pointing at low angles to the beam axis. The two-photon events have a small visible invariant 
mass Myis and a tiny transverse momentum p^is but a considerable longitudinal momentum 
along the z-axis in the common case that the two photons do not have equal energy. Due to 
these special characteristics this background can be easily reduced to a negligible level. 

The two-fermion background important for this search consists of Z/7* — > qq(7) events. 
These events tend to have a big cross-section if one or more initial state radiation photons 
(abbreviated as ISR photons) are emitted so that the effective centre-of-mass energy is 
reduced to a value near the Z-resonance (so-called radiative return events). The emission of ISR 
photons happens predominantly at small polar angles. In case of a mismeasurement or escape 
of the ISR photons through the beam pipe these events have a sizeable missing momentum 
preferentially oriented at small polar angles, close to the beam pipe. In such events the two jets 
are almost coplanar. 

The most difficult background to separate is four-fermion processes with neutrinos in the 
final state, such as W + W - ^ l^v qq and W^e^v — ► qqe^v with the charged lepton escaping 
detection. The irreducible background to this search stems from ZZ — ► vVqq (about 28 % of all 
ZZ decays) leading to a signature indistinguishable from a signal event with a Higgs mass close 
to the Z boson mass. The vector bosons are usually not produced at rest, leading to a transverse 
momentum of the two-jet system and therefore to a large acoplanarity of the jets, as in the signal 
case. Furthermore, the missing momentum vector points into the central detector more often 
than for the two-fermion case. To discriminate between this background and the signal one can 
exploit the difference in the angular distribution of the differential production cross-section. 

4.2 Preselection 

In order to reduce the amount of background data only a events fulfilling the following quality 
criteria are analysed. Prom cut No. (5) onwards, the cut values were defined using as a guide 
a simple figure of merit based on the efficiency and expected background. The following cuts 
remove almost all the two-photon background: 

(1) To reduce two-photon and accelerator induced background, track criteria are applied de- 
manding that more than 20 % of all tracks be qualified as good measured tracks [24] and 
that at least 6 of them be found. 

(2) A forward energy veto rejects events with more than 5 GeV in either the left or right 
compartment of the gamma catcher calorimeters or the silicon tungsten luminometers. 
Events with more than 2 GeV in the forward calorimeters are also removed. 

(3) The missing transverse momentum p^is snou ld exceed 1 GeV and Myis has to be larger 
than 4 GeV. 

(4) Less than 20 % of the measured visible energy -EviS should be located close to the beam 
pipe in the region | cos#| > 0.9. 

(5) The visible energy -Eyis must be less than 90 % of \fs. 

(6) It is required that the visible mass of the event should be of order Mz, i.e. 55 GeV < 
Myis < 105 GeV. An asymmetric cut around the Z mass is chosen, since with increasing 
Higgs mass Mr the Z bosons will be more and more off-shell. 

The remaining backgrounds at this stage, which are more difficult to remove, are mismeasured 
Z/7* — ► qq events, four-fermion processes with neutrinos in the final state, such as W + W" — ► 
i^vqq and W^e^v — > qqe^v with the charged lepton escaping detection (see Table 2). 
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(7) To select events that are well measured in the detector with a visible mass Myis close to 
Mi and a sizeable transverse momentum py is the following criterion is applied: Myis + 
5 x p^ ls > 

(8) A large part of the qq events and the remaining two-photon background is eliminated by 
requiring the visible transverse momentum Pyis ^ ^ GeV. 

(9) To remove backgrounds in which particles go undetected down the beam pipe, the pro- 
jection of the visible momentum along the beam axis, Pvis> 1S required to be less than 
0.294^5- 

(10) To reduce the radiative qq(7) background, the polar angle of the missing momentum vector 
must lie within the region | cos#mis| < 0.9. 

(11) The axes of both jets, reconstructed with the Durham algorithm, are required to have 
a polar angle satisfying |cos#| < 0.9 to ensure good containment. Furthermore this cut 
exploits the fact that events of the WW and ZZ background are produced according to an 
angular distribution proportional to cos 2 6. 

(12) The remaining background from Z/7* — ► qq is characterised by two jets that tend to be 
back-to-back with small acoplanarity angles, in contrast to signal events in which the jets 
are expected to have some acoplanarity angle due to the recoiling Higgs boson. Here the 
acoplanarity angle <^>acoplan is defined as 180° — <pjj where 4>jj is the angle between the 
two jets in the plane perpendicular to the beam axis. This background is suppressed by 
requiring that the jet-jet acoplanarity angle be larger than 5°. 

(13) W + W~ events with one of the W bosons decaying leptonically and the other decaying 
into hadronic jets are rejected by requiring that the events have no isolated leptons. In 
this context, leptons are low-multiplicity jets with one, two or three tracks, associated to 
electromagnetic or hadronic energy clusters, having an invariant mass of less than 2.5 GeV 
and momentum in excess of 5 GeV. In the case of a single-track candidate, the lepton is 
considered isolated if there are no additional tracks within an isolation cone of 25° half- 
angle, and if the electromagnetic energy contained between cones of 5° and 25° half-angle 
around the track does not exceed 5 % of the sum of the track energy and the electromagnetic 
energy within the 5° half angle cone. In the case of a two- or three-track candidate, 
consisting of the tracks and electromagnetic or hadronic energy clusters confined to a cone 
of 7° half-angle, the lepton is considered isolated if the sum of track and electromagnetic 
energy between the 7° half-angle cone and a 25° half-angle isolation cone does not exceed 
15 % of the lepton energy. 

For each individual centre-of-mass energy there is good agreement between the numbers of ex- 
pected background events and observed candidates after the preselection. Table 2 gives the 
number of preselected events summed over all centre-of-mass energies. Figure 1 shows the dis- 
tributions for background classes summed over all centre-of-mass energies and three arbitrarily 
scaled signal distributions (at a centre-of-mass energy of 206 GeV). The efficiencies of the pres- 
election vary on average between 39 % and 55 % for small decay widths and between 45 % and 
53 % for larger decay widths above Th = 100 GeV. 

4.3 Likelihood analysis 

To consider the changing kinematic properties of the signal hypotheses in an optimal way, five 
different likelihood-based analyses for the signal and background discrimination were applied 
after the preselection. By a likelihood analysis we denote the combination of a set of likelihood 
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input variables, a so-called likelihood, and the corresponding reference distributions of these 
variables. The reference distributions are filled with events of the specific classes for which the 
likelihood is calculated. The classes considered in this search are the two- and four-fermion 
backgrounds and the signal events. The two-photon events are negligible after the preselection. 
The search uses combinations of two likelihoods and three fixed signal mass ranges for unweighted 
reference histograms. 

To compare the kinematic properties of a selected data event to the hypothesis (Mn,Tu) 
when evaluating the likelihood, one in principle has to fill weighted signal reference distributions 
for each hypothesis (Mh,Th). This will soon lead to an unmanageable technical effort, given the 
number of hypotheses scanned. Therefore a compromise was sought in which certain kinematic 
properties of the signal were emphasised and simultaneously the number of reference histograms 
kept small. This was achieved by filling unweighted signal reference histograms. For most of the 
(Mh^h) hypotheses all signal masses were used for filling the reference histograms. This reflects 
the fact that for a very large decay width of the Higgs boson the possible values of kinematical 
variables are also smeared out over a large range. It was, however, found that the sensitivity of 
the likelihood selection (i.e. the median expected upper limits on a(Mn,Tu) x BR(H — ► -Emis)) 
could be increased further for small widths below 50 GeV by filling reference histograms with 
signal masses from 50-80 GeV and from 80-120 GeV for intermediate and heavy Higgs boson 
masses respectively. A first likelihood was designed for a signal consisting of small masses (Mr 
< 80 GeV) or large masses and a very large width (Th > 110 GeV). In this likelihood input 
variables are used exploiting the characteristics of the dominant fraction of light masses in the 
signal mass distribution. However for signal masses above Mh = 80 GeV and small or moderate 
(i.e. below 110 GeV) decay widths, the contribution of large masses dominates the signal mass 
distribution. In this case the kinematics and topology of the signal events are determined by 
the higher masses close to the kinematic limit. A second likelihood is therefore built with input 
variables optimised for such signal characteristics. In the following the choice of the inputs for 
the two optimised likelihoods are presented. 

The first three input variables are used in both likelihoods (see Figure 3) . 

(1) (l+7>(M V is ^Mz))- 1 

V(Myis = Mz) is the probability of a kinematical x 2 fit of the jet four-vectors under the 
assumption that the invariant mass of the two jets is compatible with the Z boson mass. 
The uncertainties on the measured jet energies are of the order of 5-10 GeV, while the jet 
directions are measured to approximately 1-2° [25]. This variable depends only weakly on 
the Higgs mass. For events with non-converging fit the probability is set to zero. They 
therefore accumulate at a value of 1. 

(2) -logy 32 

The Durham algorithm groups two energy flow objects i and j into a jet as long as their 
separation in phase space yij = 2 x min(Ef,Ej) x (1 — cos(0y) / Ey is is smaller then the 
cut value y cu t- The number of jets in a event is predefined to be 2, 1/32 is the value of y C ut 
where the two-jet topology of the event changes to a three-jet topology. Hence the negative 
logarithm of the so-called jet resolution parameter 1/32 is a measure for the jet topology 
being more two-jet like (large value of — log 1/32) or three-jet like (small value of — log 2/32). 

( 3 ) Pmis/Vs 

The transverse missing momentum p^is 1S one °f the most prominent characteristics of 
signal-like events, but depends very much on the Higgs boson mass. For a heavy Higgs 
boson produced close to the kinematic threshold almost at rest, the Z boson has almost no 
boost and decays into two more or less back-to-back jets. In this case the discriminating 
power of the variable is lost. 
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The next three variables (see Figure 4) complete the first likelihood, which is used for all Higgs 
masses in the domain of very large width > 110 GeV or low Higgs masses < 80GeV. 

(4) 0ACOL 

The acolinearity angle <^acol of the two-jet system is obtained by subtracting the three- 
dimensional angle between the reconstructed jet-axes from 180°. Events containing a low- 
mass Higgs boson exhibit on average a larger acolinearity than the background. 

(5) | cos 6* | 

The Gottfried- Jackson angle 6* , is defined as the angle between the flight direction of the 
Z boson in the laboratory frame and the direction of the decay products of the Z boson 
boosted into the Z boson rest-frame. The variable tends to have smaller values for the 
signal. 

(6) -log y 2 i 

The variable — log y2i is analogous to — log 1/32 and measures the compatibility of the event 
with a two-jet topology. Two-jet events tend to accumulate at small values of — log 2/21 • 

The last three variables (see Figure 5) tune the second likelihood to become more sensitive for 
large Higgs boson masses and small to moderate widths. 

(7) £j M E a T 7v^ 

The variable measures the energy of the most energetic of the two jets. This is on 

average higher for the four-fermion background, due to the boost of the W and Z pairs, 
whereas heavy Higgs bosons and a Z boson are produced at rest. 

(8) Rp ti 

This variable is the significance of the acoplanarity between the two jets, taking into account 
detector resolution and acceptance. The discrimination power is enhanced by weighting 
the acoplanarity with the average jet polar angle, since transverse jet directions are more 
precisely measured at large polar angles. Signal events tend to have a more significant 
acoplanarity and thus larger values of Rp ti than background. The precise definition of Rp ti 
can be found in the OPAL analyses of ZZ — > qquV events [26]. 

(9) (Myis + M M is)/(M V is - Mmis) 

This variable, described in [27], uses two strongly correlated quantities, the invariant miss- 
ing mass Mmis and the visible mass of the event Myis- Depending on the mass recon- 
struction accuracy it can have positive or negative values. The signal distribution of this 
variable is broader and accumulates at higher values than for the two- and four-fermion 
events, which are distributed more narrowly around the origin. 

From the two likelihoods and three ranges of signal masses filled in the reference histograms 
one has six analyses to search for the different hypotheses in Mh and Th- The study of the 
median expected a(Mn,Tn) x BR(H — ► E'mis) shows that five of these six are sufficient to have 
an optimally efficient analysis for each signal hypothesis characterised by Mh and Th (see Figure 
2) in the range studied. Likelihood 1 was not used with the reference distribution filled for the 
signal mass range of 80 to 120 GeV. Figure 6 a) to c) and g) to h) display examples for the 
likelihood distributions of all five analyses used. In the histograms the events selected at all 
five centre-of-mass energies are added up, although each centre-of-mass energy was evaluated 
separately in the limit setting as explained in Section 5.1. The appropriate likelihood was cal- 
culated for each background, data and signal event. In case of a signal event it was added to 
the histogram with the weight defined in Equation 7. The number of expected signal events is 
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normalised according to Equation 6. The use of different analyses gives rise to varying shapes of 
the likelihood distributions of the background. Also the various shapes of the signal likelihood 
distributions are visible for different Mh and Th- Since the form of the likelihood distributions 
for signal and background can yield additional information in the limit calculation, only a loose 
cut is applied in the likelihood selection, requiring a signal likelihood larger than 0.2. 

4.4 Correction of background and signal efficiencies 

A correction is applied to the number of expected background events and the signal efficiencies 
due to noise in the detectors in the forward region which is not modelled by the Monte Carlo. 
The forward energy veto used in the preselection can accidentally be triggered by machine back- 
grounds. The correction factor is derived from the study of random beam crossings, and applied 
individually for each year of data taking. Random beam crossing events were recorded when no 
physics trigger was active. The fraction of events that fail the veto on activity in the forward 
region is below 3.4% for all runs analysed. The detailed breakdown of the fraction of accidentally 
vetoed events is given in the last column of Table 1. 

4.5 Systematic uncertainties 

A possible signal in the data would reveal itself by altering the shapes of the distributions of 
the discriminating variables. Thus a systematic deviation in the description of a reconstructed 
observable between Standard Model Monte Carlo and a data sample in which the signal is absent, 
could wrongly be attributed to the presence of a signal. 

The systematic uncertainties in the Monte Carlo description of the kinematic event variables 
are studied in two control samples at a centre-of-mass energy of 206 GeV. In the first control 
sample, called two-fermion control-sample in the following, radiative returns contributing to the 
qq(7) processes with photons detected at large angles are selected and the tagged ISR photon is 
removed from the event in Monte Carlo and data. This creates a qq-like topology with missing 
momentum at large angles. The second control sample, called four-fermion control-sample, is 
obtained by selecting W + W~ — > qqlv events and removing the identified isolated lepton from 
the events in Monte Carlo and data. After this procedure these two control samples possess a 
topology very similar to signal events. For all kinematic variables x of the preselection and the 
likelihood selection the mean x and the width of the distribution (RMS) are compared between 
Monte Carlo and data, for the two-fermion and four-fermion control-samples. The observables 
in the two-fermion and four-fermion Monte Carlo are then modified separately according to 
x mc W = ( x mcP ~^mc) x R rmsmc +^ data - Then all five likelihood selections are repeated 
separately and the relative change in the number of selected events compared to the unmodified 
case is taken as the systematic uncertainty. 

Since the analysis labelled Al in Table 3 is used over a large range of the search plane 
(see Figure 2), the systematic uncertainties on the background determined in this analysis Al 
are taken as an estimate for the background for all analyses. To determine the effect of the 
systematic uncertainty on the signal efficiencies, one has to take into account the fact that 
the kinematic properties of the signal depend on the assumed Higgs mass and decay width. 
Twelve representative hypotheses are studied with Mh chosen to be 20, 60 or 110 GeV and 
Th taking values of 5, 20, 70 and 200 GeV. For these hypotheses the signal Monte Carlo is 
modified according to the four-fermion correction factors, representing the dominant remaining 
background after the cut on the signal likelihood. The relative change in selected event weights 
compared to the unmodified case is then taken as an estimate of the systematic uncertainty on 
the signal efficiency for a given hypothesis. The root-mean-square of all twelve hypotheses is 
applied as an (MH,rn)-independent estimate for the whole search area and for all centre-of-mass 
energies (see Table 4). 
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The W pairs are very effectively reduced in the preselection by the isolated lepton veto. Due 
to the importance of this veto the uncertainty from the lepton isolation angle and the vetoed cone 
energy is studied in the following way. The half-cone angle of the outer cone is increased and 
decreased by two degrees, following the studies in [28], and the relative effect on the selection 
determined. Furthermore the cone energy is varied by 7.4 % and the analyses are repeated. 
The value of the cone energy rescaling is determined by the relative deviation of the mean of the 
measured energy of the lepton candidates in the inner cone between data and Monte Carlo in the 
W + W _ — > qqlu sample. For signal efficiencies an analogous study was performed at the twelve 
points described above. Both results for the relative change of the selection for the cone opening 
half-angle and cone energy variation are added in quadrature and the root-mean-square of the 12 
(Mh,Fh) hypotheses was taken to yield the total uncertainty associated with the isolated lepton 
veto (see Table 4). 

The theoretical prediction on the cross-section for the two- and four-fermion processes adds 
an uncertainty of 2% to the background uncertainty [29]. Finally, the uncertainty due to the 
limited Monte Carlo statistics is evaluated. 

Table 5 summarises the results of the studies. All uncertainties are assumed to be uncorre- 
cted and the individual contributions are added in quadrature to obtain the total systematic 
uncertainties on the background expectation and signal efficiency The dominant systematic 
uncertainties on the signal efficiency arise from the description of the kinematic variables. The 
background expectation is more affected by the uncertainty in the isolated lepton veto, as the 
main contribution of the background stems from four-fermion processes. But the uncertainty 
associated with the description of the kinematic variables is of similar magnitude. The limits 
quoted in Section 5.1 were calculated including the uncertainties of Table 5. To estimate the 
extent to which the limits depend on the size of the systematic uncertainties, the limit calculation 
was repeated doubling the systematic uncertainties. A comparison of the limits with single and 
double systematic uncertainties, done at similar representative points as used for the systematic 
studies, showed that the excluded cross-sections typically decrease between a half and one and 
a half percent. A maximal reduction of 2.1 % was found. 

5 Results 

The results of the search using each of the five different likelihood selections, labelled A1-A5, after 
a cut on the likelihood larger than 0.2 are summarised in Table 6, which compares the numbers of 
observed candidates with the Standard Model background expectations. The data are compatible 
with the Standard Model background expectations. The remaining four-fermion background 
consists predominantly of W pairs, representing roughly three quarters of the background at 
energies above the Z pair threshold. Figure 7 shows examples for signal efficiencies. For small 
decay widths the dependence on the centre-of-mass energy is weak up to Mh ~ 80 GeV, and 
for large widths it is weak up to the kinematic limit. Because of the centre-of-mass energy 
dependence of the Bjorken cross-section of a Higgs boson with mass Mh, the lower centre-of- 
mass energies contribute more significantly to the sensitivity for lighter Higgs bosons (see e.g. 
in [8]). For very light Higgs bosons the efficiency is moderately reduced by the preselection cuts 
demanding a sizeable amount of missing energy. In the case of broader Higgs resonances with 
high mass, one observes a generally enhanced efficiency since the chance of selecting events from 
the low mass tail compensates the suppression due to the falling production cross-section of a 
heavy Higgs boson. 
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5.1 The upper limits on the production cross-section times branching ratio 

Upper limits are calculated on the model-independent cross-section a(Mn, Th) x BR(H — > Emis) 
scaled to \fs = 206 GeV. As the likelihood distributions are only loosely cut, one can use not 
only the information from the integral number of selected events (Table 6) but also from the 
shape in a likelihood ratio [30] to set more sensitive upper limits. For each centre-of-mass energy 
separately, each bin with a signal likelihood larger than 0.2 in the distributions of expected 
signal, background and selected data is treated as a search channel. For each centre-of-mass 
energy the number of expected signal events is scaled to the total cross-section (Equation 6). As 
with the analysis described in [31], the likelihood distributions are given as discriminating input 
to a limit program [32]. A likelihood ratio is used to determine the signal confidence level, CLs, 
defined in [30, 32], which excludes the presence of a possible signal according to the modified 
frequentist approach [32]. Additionally the program calculates the median upper number of 
signal events that could be excluded at 95% confidence level (CL). This number is then scaled 
to the total cross-section at the centre-of-mass energy of 206 GeV for each (Mh^h) hypothesis. 
The systematic uncertainties on the background expectations and signal selection efficiencies are 
included according a generalisation of the method described in [33]. 

A very fine scan of the (Mn,rH)-plane was performed by simulating the spectra of Higgs 
bosons with a mass Mh from 1 to 120 GeV and widths Th starting at 1 GeV up to 3 TeV. The 
Higgs boson mass was simulated in steps of <5Mh = 1 GeV. Simulated values of Th are spaced 
in steps of 1 GeV up to 5 GeV. A spacing of 5Tn = 5 GeV is chosen from Th = 5 GeV to Th = 
750 GeV. Above this value steps of 5Tn = 50 GeV are adopted up to the maximal Th of 3 TeV. 

Examples of the projections of the observed upper cross-section limits together with the 
median expected upper limits and the corresponding one and two standard deviation bands 
on the expected limits are displayed in Figure 8 for some choices of Th- Above a width of 
300 GeV the exclusion plots look quite similar to the example displayed in Figure 8i) because 
the excluded limits do not change very much. The observed limits for Th > 60 GeV are well 
contained in the one standard deviation bands on the expectations and generally do not exceed 
two standard deviations except in Figure 8 a) at Mh = 114 GeV. The discontinuities in the 
graphs correspond to changes in the analyses. As one can observe, below Th < 40 GeV the 
analysis are changed more often. Therefore the chance is higher that in a few bins there are 
statistical fluctuations in the selected data, that lead to a deviation of more than one standard 
deviation around the median. Also the data selected are highly correlated, as one can see for 
example in the upward fluctuation around Mh = 114 GeV visible in Figure 8 a)-c). All results 
for the observed upper limits on ct(Mh,Fh) x BR(H — > Emis) are summarised in a contour plot 
(Figure 9) in the scanned (MH,rn)-plane. Above Th = 200 GeV the observed upper limits are 
in the range of 0.15 pb to 0.18 pb for all Mh and vary very little. For such large Th the recoil 
mass distribution of the Higgs tends to be more and more uniformly stretched out over the mass 
range explored. There is not much difference in the selection of signal events for a Higgs boson 
with e.g. a width of 400 GeV or 600 GeV in the considered range of Higgs masses. This prevents 
any specific discriminating kinematical properties from being assigned to the expected signal as 
signal masses of a broad kinematical range are selected with roughly equal probability. Therefore 
only one likelihood analysis is used in this part of the search area, selecting the same subset of 
data and background. Since the upper limit on the model-independent cross-section refers to a 
production cross-section at a centre-of-mass energy of 206 GeV, it must become independent of 
(Mh,Fh) for an extremely large Th- In this case the shape of the Higgs signal would just be 
a box, weighted with the production cross-section from 1 GeV to the kinematic limit of about 
115 GeV. The data are then compared to an approximately constant signal expectation. Hence 
the upper limit on the cross-section is approximately independent of the (Mn,Tn) hypothesis at 
a value of roughly 0.16 pb. For resonances with a decay width smaller than 200 GeV there are 
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regions where the limits are below 0.15 pb or even 0.1 pb for Mh between 60 and 74GeV. In this 
mass range the number of data events selected is smaller than expected. Above Mh of 85 GeV 
the upper limits become larger than 0.2 pb and rise considerably for small widths below 40 GeV 
(see Figure 8 a-e). This is due to the fact that the Higgs mass approaches the kinematic limit 
and the likelihoods which rely on kinematical variables like p^is l° se discrimination power. A 
maximal value of 0.57pb is observed for Mh of 114 GeV and Th of 1 GeV corresponding to a 
circa two-cr excess in the data. 

It should be kept in mind that no optimisation of the search has been performed for Th 
below 5 GeV. In the region of heavy Higgs boson mass > 105 GeV and small width a search 
using recoil mass spectra would be more sensitive. Therefore this region is more sensitively 
covered by searches that have been performed by the LEP experiments documented in [2]. 

5.2 Interpretation of the result in the stealthy Higgs scenario 

Interpreting the width Th of a Higgs boson according to Equation 5, and setting m p hi n to zero, 
it was possible to set limits on u in the stealthy Higgs scenario. A range from u> = 0.04 to 
ijj = 24.45 was probed. The excluded regions are shown in Figure 10 at 95 % confidence level 
(CL) in the to-Mn parameter space. To illustrate the Higgs boson width according to Equation 5, 
contours of fixed Th corresponding to a given mass Mh and coupling lo are added to the plot. 
The maximum excluded invisible width is about Th = 400 GeV for Higgs boson masses < 35 GeV, 
decreasing slowly to Th = 115 GeV for Mh = 100 GeV. The minimal exclusion of u = 0.04 is 
observed at Mh = 1 GeV and the maximal exclusion is to = 5.9 for Mh = 73 GeV. For u> between 
0.04 and 0.59 a Higgs mass from 1 to 103 GeV could be excluded. The maximal excluded Higgs 
mass was 103 GeV for width between 1 and 3 GeV, compared with the expected exclusion of 
106 GeV. 

The results presented in this study extend the previous decay-mode independent searches 
for new scalar bosons with the OPAL detector [5] to regions of larger couplings and higher 
Higgs boson masses. In [5] an interpretation within the stealthy Higgs model yielded a maximal 
excluded coupling uj for masses around 30 GeV, where lo was excluded up to lo = 2.7. That study 
excluded Higgs boson masses up to Mh = 81 GeV. It should be pointed out that the decay-mode 
independent searches also studied Higgs widths between 0.1 and 1 GeV and therefore cover the 
gap between searches within scenarios assuming a narrow decay width of the invisibly decaying 
Higgs boson [2] and the search presented in this paper up to Mh = 81 GeV. 

6 Conclusions 

A dedicated search was performed in the channel e + e~ — > HZ with Z — > qq and the non- 
standard Model decay H — ► -Emis final state allowing for invisible decay widths of the Higgs 
boson from 1 GeV up to 3TeV. The data taken by the OPAL detector at LEP above the W 
pair threshold were analysed. No indication for a signal was found and upper limits were set on 
£7(M H ,r H ) x BR(H -»■ £mis)- The maximal upper limit is 0.57 pb at M H = 114GeV and T H = 
1 GeV. Over the scanned region of the (MH,rn)-plane upper limits are generally of the order of 
0.15 pb, especially for large values of Th > 400 GeV or Higgs boson masses < 85 GeV. 

The limits were interpreted in the stealthy Higgs scenario assuming the presence of a large 
number of massless singlet states. Limits were calculated on the coupling lo to a hidden scalar 
sector of the Higgs boson with a given mass Mh. A large part of the parameter plane kinemati- 
cally accessible with LEP 2 was excluded extending a previous exclusion published in [5] . Values 
for lo between 0.04 (Mh = 1 GeV) and 5.9 (Mh = 73 GeV) were excluded, and for certain values 
of lo Higgs boson masses are excluded up to Mh = 103 GeV. The possible non-detection of a 
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light Higgs boson at the LEP searches due to non-Standard Model invisible Higgs boson decays 
is therefore restricted to the case of extremely large decay widths > 400 GeV. 
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binned y/s 


nominal y/s (GeV) 


year 


int. luminosity (pb 1 ) 


accid. veto (%) 


> 180-186 


183 


1997 


40.0 


3.37 


> 186-193 


189 


1998 


199.8 


2.24 


> 193-198 


196 


1999 


70.4 


2.53 


> 198-203 


200 


1999 


112.0 


2.96 


> 203-209 


206 


2000 


206.9 


2.22 



Table 1: Breakdown of the analysed integrated data luminosities according to the centre-of-mass 
energies. The data was binned in five nominal centre-of-mass energies. The last column states 
the reduction of the signal efficiencies and expected background rates due to accidental triggering 
of the forward energy veto in the preselection, which is not modelled in the Monte Carlo. 



cut 


77 


qq( 7) 


4-fermion 


total SM 


data 


(l)-(5) 


48795 


15639 


4880 


69314 


74178 


(6) 


148 


10359 


1394 


11901 


11779 


(7) 


62 


9128 


1336 


10526 


10472 


(8) 


44 


4897 


1167 


6108 


6264 


(9) 


33 


1061 


964 


2058 


2116 


(10) 


18 


425 


895 


1338 


1387 


(11) 


18 


423 


879 


1320 


1368 


(12) 


4 


68 


820 


892 


899 


(13) 


4 


60 


441 


505 


498 



Table 2: Expected number of Standard Model background events after the preselection nor- 
malised to a data luminosity of 629.1 pb -1 . The total SM background after preselection is 
expected to be 505 ± 5(stat) ± 21(syst). The contributions of the different subclasses are bro- 
ken down in column two to four for the two-photon, two-fermion and four-fermion processes 
respectively. 



label 


likelihood 


analysis 

reference mass range (GeV) 


background uncertainty 
kinematic var. isol. lepton veto 


Al 


1 


1-120 


2.4% 


2.4% 


A2 


2 


1-120 


1.6% 


2.3% 


A3 


1 


50-80 


1.0% 


2.5% 


A4 


2 


50-80 


1.6% 


2.6% 


A5 


2 


80-120 


1.1% 


1.5% 


choice for uncertainty 


2.4% 


2.4% 



Table 3: Results of the study of systematic uncertainties of the expected background for the 
five kinds of analyses, labelled A1-A5, used in the search (see Figure 2) at a centre-of-mass 
energy of 206 GeV. Since the analysis labelled Al covers the largest part of the search area, its 
uncertainty was chosen as representative uncertainty on the background due to the uncertainty 
in the kinematic variables and the isolated lepton veto at all centre-of-mass energies. 
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signal hypothesis 
M H (GeV) T H (GeV) 


efficiency uncertainty 
kinematic var. isol. lepton veto 


zO 


5 


0.6 7o 


0.6 7o 


on 


on 


U.I A) 


U. 1 A) 


20 


70 


n 3 % 

u.o /u 


7% 

U. ( /u 


20 


200 


0.1% 


0.7% 


60 


5 


0.7% 


0.8% 


60 


20 


0.7% 


0.8% 


60 


70 


0.2% 


0.8% 


60 


200 


0.3% 


0.7% 


110 


5 


5.5% 


0.7% 


110 


20 


2.9% 


0.8% 


110 


70 


1.3% 


0.8% 


110 


200 


0.1% 


0.8% 


all M u 


and Th 


1.9% 


0.7% 



Table 4: Results of the study of systematic uncertainties in twelve representative (Mu,Tu)-points 
at a centre-of-mass energy of 206 GeV. For each source the root-mean-square of the individual 
uncertainties in the twelve points was taken to get an (Mn,Tn) independent estimate of the 
uncertainty at all centre-of-mass energies. 



source 


background uncertainty 


efficiency uncertainty 


kinematic variables 


2.4% 


1.9% 


isolated lepton veto 


2.4% 


0.7% 


limited MC statistics 


1.0% 


0.2% 


prediction 2- and 4-f cross-sect. 


2.0% 




total uncertainty 


4.1% 


2.0% 



Table 5: Results of the study of systematic uncertainties of the background for the five analyses 
(see Table 3) and of the signal efficiencies in twelve representative (Mn,Tu)-points (see Table 4) 
at a centre-of-mass energy of 206 GeV. The total uncertainty on background expectation and 
signal efficiency is applied at all centre-of-mass energies and for all (Mn,Tu) hypotheses. 



label 


reference masses 


likelihood 


2-fermion 


4-fermion 


total SM 


data 


Al 


1-120 GeV 


type 1 


11 


374 


385 ± 4 ± 16 


369 


A2 


1-120 GeV 


type 2 


3 


378 


381 ± 4 ± 16 


370 


A3 


50-80 GeV 


type 1 


5 


315 


320 ± 3 ± 13 


305 


A4 


50-80 GeV 


type 2 


2 


315 


317 ± 3 ± 13 


310 


A5 


80-120 GeV 


type 2 


8 


247 


255 ± 3 ± 11 


253 



Table 6: The likelihood selection of events with a signal likelihood exceeding 0.2 according 
to the different search strategies. The individual contributions to the total Standard Model 
background of two-fermion and four-fermion background is broken down in the second and third 
column respectively. For the total Standard Model background the statistical and the systematic 
uncertainty is also given. 
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Figure 1: Distribution of the preselection variables after the preselection cuts (l)-(5). All classes 
of Standard Model background and data are added for all analysed centre-of-mass energies. The 
distributions of three arbitrarily scaled signal hypotheses at = 206 GeV are displayed as open 
histograms. 
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Figure 2: In total five analyses were used to cover the plane of hypothetical Higgs mass and decay 
width pairings. The analyses differ in whether the first or second likelihood was used (denoted 
by the number in the cell) and what signal masses where used in filling the reference histograms 
(depicted by the shading of the cell). The pattern resulted from an optimisation starting with 
= 5 GeV up to 50 GeV. Below Th = 5 GeV the pattern was simply continued and not 
optimised anymore. Above 50 GeV a simple continuation of the pattern was found and proved 
to be sufficiently sensitive. 
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Variables used in likelihoods 1 and 2 

• OPAL 

I I 2-fermion 
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-■ (M H ,r H )=( 80,15) (GeV) 
— (M H ,r H )=( 40,20) (GeV) 




(1+P(M VIS =M 7 )) 




Figure 3: Distributions of the likelihood variables. All classes of Standard Model background 
and data are added for all centre-of-mass energies analysed. The distributions of three arbitrary 
scaled signal examples at yfs = 206 GeV are displayed as open histograms. The variables shown 
contribute to likelihood 1 and 2 as they exploit general properties of the signal signature. 
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Figure 4: Distributions of the likelihood variables. All classes of Standard Model background 
and data are added for all centre-of-mass energies analysed. The distributions of three arbitrary 
scaled signal examples at y/s = 206 GeV are displayed as open histograms. The variables shown 
are combined with the ones of Figure 3 to construct the likelihood 1 used in a general search at 
different Mh and Th- 



24 



Variables used in likelihood 2 

• OPAL 
I I 2-fermion 
I I 4-fermion 

(M H ,r H )=(110,10) (GeV) 
-- (M H ,r H )=( 80,15) (GeV) 
— (M H ,r H )=( 40,20) (GeV) 
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(M VIS+ M MIS )/(M VIS -M MIS ) 



Figure 5: Distributions of the likelihood variables. All classes of Standard Model background 
and data are added for all centre-of-mass energies analysed. The distributions of three arbitrary 
scaled signal examples at yfs = 206 GeV are displayed as open histograms. The variables shown 
have a larger discrimination power for a heavier Higgs boson and contribute with the variables 
of Figure 3 to the second likelihood. 
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Figure 6: Examples of some of the likelihood selections. Figure 6 e), f), h) corresponds to analysis 
Al (as labelled in Table 3 and 6), a) to A3, b) to A4 and c), d) to A5. The OPAL data and the 
expected 2-fermion and 4-fermion background are added for all analysed centre-of-mass energies. 
The signal hypothesis in the hatched histograms is normalised to the number of expected signal 
events and added to the background. The open histograms display the shapes of scaled signal 
distributions. 
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Figure 7: Examples for the selection efficiency after a cut on the signal likelihood greater than 
0.2 versus the Higgs mass Mh as function of the assumed decay width Th at the different \/s. 
The error is the binomial error on the selected event weights and smaller than the markers. 
Lines are added to guide the eye. A signal in the range of 80 to 90 GeV suffers from a drop in 
the efficiency due to the relatively large remaining W- and Z-pair backgrounds. For a smaller 
widths Th the efficiency to detect a relatively heavy (above 100 GeV) and more Standard Model 
like Higgs boson is more restricted by the available \/s. For a large Th signal hypothesis, the 
kinematic distributions of events and the distribution of weights assigned to these events are 
broader. Therefore it is more likely to select a larger fraction of the event weights leading to 
more uniform efficiency, that does not depend very much on the centre-of-mass energy. 
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Figure 8: The model independent upper limits at 95 % CL on the production cross-section times 
branching ratio, <7 pro d X BR mv , scaled to a centre-of-mass energy of 206 GeV for Higgs mass Mr 
and some examples of the Higgs decay width Th. The discontinuities in the limits reflect the 
changes in the analysis used at this mass (see Figure 2). 
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Figure 9: The exclusion contours at 95 % CL on the model independent production cross-section 
times branching ratio, a pro d x BR mv , for Higgs boson mass Mh and the Higgs boson decay width 
Th up to 3 TeV (note a change in logarithmic scale below Th = 5 GeV for better visibility). Solid 
and dashed lines delimit areas of excluded upper limits. Cross-sections times branching ratio 
between 0.07 pb and 0.57 pb are excluded with the OPAL data above yfs = 183 GeV. 
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Figure 10: The exclusion contours at 95% CL on the Higgs mass Mh and the Higgs-phion 
coupling u of the stealthy Higgs model. The values of u are related to the decay width Th via 
Th(Mh) = Fsm(-^h) + 32*7]- Mh ' * n ^ e case °^ mass ^ ess phions (see Equation 5). Contours of 
fixed Th are also shown in the plot as dashed lines. 
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